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NOTATION 


R 
A, Expanded blade area, Ue ear 
AL/Ay Propeller expanded area ratio, (22/n)s) (c/D) dx 
h 
(Ac/Ay),. Keller's minimum expanded area ratio for eliminating back 
bubble cavitation, (2.6+0.6Z)K1/{op glJ2+(0.7n)?]}+K 
Ay Disc area, 1D2/4 
Ap Estimated propeller projected area, 
[1.067-0.229 (P/D) JA, 
a(x) Area of section, 2c (x)t (x) spt (xx, lox 
B(x) Distance of CG from hub face, 


y cos 4 +X sin($-8,) x R tan($-8,) x4D,,/2 


(c/R)) 6s (C/R) o¢ Chord lengths measured from leading edge and trailing edge 
of blade to propeller reference line 


Cy Section drag coefficient 
Cro Frictional resistance of section ? 
a6 Center of gravity 
C Blade section lift coefficient 
Cp Power loading coefficient, Pyflo/2)=R2V,3} 
Cos Power loading coefficient based on ship speed, 
Py/L(e/2)aR-V9); calculated im (1 +e/tanB, ) 
(dCpcq/dx) ax : 
oa Thrust loading coefficient, T/L(o/2}#RV,"] 
Crs Thrust loading coefficient based on ship speed, 


T/L (p/2)aR2V2); calculated _ (1 - etanB,) 
(dCy. ;/dx) dx h 


CThp Power loading coefficient, i (1 - W) (1 - etang,) 
(dCy<;/dx) dx h 
i 
c Propeller blade chord length, c(x) 
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Inviscid power loading coefficient, 
(42/2, )xGL (I=, #4, /2V] 


Inviscid thrust loadine coefficient, 
426[(x/d, )-U,/2V] 


Propeller diameter 
Hub diameter 


Parameter for calculating the fluctuating 
angles of attack, 1/[1+2ntan(8,-8)/C, ] 


Camber 
Acceleration due to gravity 


Nondimensional circulation about a blade 
section 1r/(2nRV) 


Spacing between fillets 
Spacing between blades at the hub 
Static head at propeller shaft centerline 


Moment of inertia of blade section about 
x and y axes 


Advance coefficient, V(1-W7)/(nd)=V,/(nD) 
Ship speed advance coefficient, V/(nD) 


Kellers' constant for predicting minimum 
blade area of propeller (see p. 22) 


25) 
Torque coefficient, Q/(m D 

204 
Turust coefficient, T/(pn OD ) 


Propeller lift distribution per unit span 
for finite element stress calculations 


Moment of blades, see page 34 
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Moment due to thrust and torque 


Moment parallel and perpendicular to 
the nose - tail line 


Propeller revolution per unit time 


Estimated prcreller pitch ratio at 0.7 
radius, 0.7ntanB, in program 


Delivered power at propeller, 27Qn 

Effective power 

Shaft power 

Propeller torque 

Propeller tip radius 

Propeller revolutions per minute 

Propeller local radius 

Propeller hub radius =. 
Local position along the section chord - 
Propeller thrust 


Propetler blade maximum thickness t(x), 
thrust deduction fraction 


Chordwise distribution of section thick- 
ness (NACA 66 modified thickness form is 
used) 

Axial induced velocity at lifting line 
Tangential induced velocity at lifting line 


Ship speed 
Speed of advance of the propeller, V(1-We) 


Local velocity along the x axis at any field 
point 


Inflow velocity at each propeller section, 
V =, )+ al VJ4+ dom z/e < 


Axial velocity from sources other than the 
propeller wake (1-w,) 


Weight of blades 
Weight of hub 
Propeller weight 


Circumferential mean wake fraction at each 
radius calculated from wake survey 


Tangential velocity from sources other than 
the propeller wake (1-w,) 


Propeller effective wake fraction as determined 
from thrust identity from self propulsion 
experiment 

Volume mean wake fraction 

Propeller wake fraction, 1-[(1-we)/(1-w (1-4, 
Nondimensional radial distance, r/R 
Nondimensional hub radius, r/R 

Nondimensional distance along section chord, r/c 
Number of blades 


Propeller rake 


Total rake, rake plus skew induced rake 
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99,7 


Section ideal angle of attack, 1.54¢C, 
for NACA a=0.8 meanline in two 
dimensional flow 


Maximum fluctuating angle of attack, 
a,~(-08)F(x) 


Minimum fluctuating angle cf attack, 
a, -(+48)F (x) 


Advance anqle of a propeller blade section 


Hydrodynamic flow angle of a propeller blade 
section 


Circulation about a propeller blade section, 


Section drag-lift ratio, C)/e, 


Propulsive efficiency, Pe/Pya(1-t)Cr/Cp¢ 


Estimated propeller efficiency,Co3/Cy¢ 
Blade rake angle in deqrees 
Blade skew angle in degrees 


Advance ratio of propeller based on ship speed, 
V/(mnD) 


Water density 
Density of propeller material 
Pitch angle 


Section cavitation number, 2gH/V* 


Burrill cavitation number, : 
2gH/[{V(1-w_9 7)} 2 +{0.7nnD} J] 


Burrill thrust loading coefficient, 
T/L(V(1-W 2g. 7)} 2 #{0.7nnD) 2 J 


ABSTRACT 


This report presents a computer program that can be used for 
the preliminary design and to predict the performance of single 
screw propellers when designed for a prescribed hydrodynamic 
pitch distribution. This design program is based on lifting line 
theory as developed by Lerbs for moderately-loaded finite-bladed 
propellers. Stress calculations using beam theory, propeller 
weight, moments of inertia, and center of gravity for specified hubs 
are also made for each design, which take into account the effect 
of blade skew and rake. Design calculations can be made for given 
design speeds where the corresponding shaft power is computed 
(thrust option), or for the case when the shaft power is specified 
and the corresponding speed is computed (power option). A FORTRAN 
listing of the program, developed to run on the computers at the 
David W. Taylor Naval Ship Research and Development Center (DTNSRODC) 
is presented as well as the input and cutput obtained for two sample 


designs, one using the thrust option, and the other the power option. 
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ADMINISTRATIVE INFORMATION 
This work was sponsored by the Naval Ship Systems Command, 
SHIPS 034, and carried out under the David W. Taylor Naval Ship 
Research and Development Center (DTNSRDC) Work Unit No. 1524-462, 
Task 15942. 


INTRODUCTION 

The David W. Taylor Naval Ship Research and Development 
Center (DTNSRDC), Carderock Laboratory, was requested by the 
Naval Ship Systems Command (NAVSHIPS) to develop a computer 
program that can be used as an aid in determining optimum 
values of propeller design parameters such as diameter, rpm, 
and blade number for naval vessels. This report has the 
objective of presenting a computer program that can be easily 
used for the preliminary design and for prediction of the 
performance of propellers applicable to specific snips. 

The main portion of the computer program makes preliminary 
propeller design calculations and performance predictions 
using Lerbs' lifting line moderately-loaded finite-bladed 
propeller theory of References 1 and 2. Once the lifting 
1. Lerbs, H.W., "Moderately Loaded Propellers with a Finite 

Number of Blades and an Arbitrary Distribution of Circulation," 

Transactions of the Society of Naval Architects and Marine 

Engineers, Vol. 60, p 73-117, 1952 
2. Morgan, W.B. and Wrench, J.W., Jr., "Some Computational 


Aspects cf Propeller Design," Methods in Computational Physics, 
Vol. 4, Academic Press Inc., New York, p 301-331, 1965 


nS 


oil 


Tine calculations for the propelier have been obtained, a 
computer program based on lifting surface theory must be 

used to determine the final pitch and camber of the propeller’. 
Numerous experimental checks’ on the design procedure used 
have verified this procedure. 

A computer program for the preliminary design of propellers 
having a prescribed pitch or circulation distribution was 
previously published in Reference 5. One of the most 
important new features of the new computer program is the 
power option which computes the speed automatically when the 
design power is specified. Other features included in the 
present computer program follows: 

1. Only one set of basic input data is required for the 
preliminary design of a series of propellers with rpm, 
blade number and expanded area ratio varied. Any 
number of basic sets of data can be specified in 


the computer program. 


ee 


3. Morgan W.B., Silovic, Vladimir and Denny, Stephen 8., "Propeller 


Lifting-Surface Corrections," Transactions of the Society of Naval 
Architects and darine Engineers, Vol. 76, p 309-347, 1968 


4. Boswell, R.J., "Design, Cavitation Performance and Open-Water 


Performance of a Series of Research Skewed Propellers," Naval 
Ship Research and Develpment Center Report 3339, March 1971 


5. Haskins, E.W., "Calculations of Design Data for Moderately Loaded 


Marine Propellers by Means of Induction Factors," Naval Ship 
Research and Development Center Report 2380, September 1967 
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2. The propeller stress, based on beam theory®, modified 
to account for the effect of rake and skew, is computed for 
each design. An input option is used to make stress calcula- 
tions for a linear or nonlinear distribution of blade skew. 

2. Parameters required to make blade surface cavitation 
checks using the Burrii] cavitation charts of Reference 7 and 
Brockett's incipient cavitation diagrams of Reference 8 are 
calculated on the computer. In addition, Keller's method? of 
predicting the minimum expanded area ratio based on back bubble 
type cavitation, is calculated. The estimated propeller 
weight including the hub, and input parameters required for 
the lifting surface computer program used to calculate the 
final pitch and camber of propeller designs are some of the 


other parameters calculated. 


6. Eckhart, M.K. and Morgan, W.B., "A Propeller Design Method," 
Transactions of the Society of Naval Architects and Marine 
Engineers, Vol. 63, p 305-370, 1955 

7. Burrill, L.C. and Emerson, A., "Propeller Cavitation: Further 
Tests on 16 Inch Propeller Models in the Kings College Cavita- 
tion Tunnel," Transactions of the North East Coast Institution 
of Engineers and Ship Builders, Vol. 78, p 295-320, 1963-64 

8. Brockett, Terry, "Minimum Pressure Envelopes for Modified 
NACA-66 Sections with NACA a=0.8 Camber and BUSHIPS Type I 
ae II, Sections," David Taylor Model Basin Report 1780, 

9. Keller, J. Auf'm, "Enige aspectem bij het ontwerpen van 
Scheepsschroeven," Schip en Werk, Ho. 24, p 658-662, 1966 
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PROPELLER LIFTING LINE THEORY 

The availability of high speed digital computers has made 
it now possible to make use of more adequate mathematical models 
to represent the hydrodynamic action of marine propellers. 

The computer has also released the designer from performing 
laborious computations involvad in present design methods. 

The lifting line theory used in the computer program 
presented here is the theory developed by Lerbs (Reference 1) 
for moderately-loaded finite-bladed propellers. This lifting 
line design method discussed by Cox and Morgan in Reference 10 
also permits the designer to accurately account for propeller 
parameters such as number of blades, hub size, radial blade- 
loading, and wake distribution. Since this theory is discussed 
in detail in References 1 and 2, only some of the assumptions 
made in developing the theory are mentioned here. The theory 
developed by Lerbs' considers the influence of the induced 
velocities on the shape of the helical vortex sheet at the 
lifting line, but neglects the effect of centrifuga! forces 


and the contraction of the slip-stream. In addition, the 


change in shape of the vortex lines is neglected in the axial 
direction, indicating that each vortex line is assumed to be of 
constant pitch. However, the vortex sheets are not necessarily 
helical surfaces since the pitch may vary along the radius. 


10. Cox, G.G. and Morgan, W.B., "The Use of Theory in Propeller 
Design," Marine Technology, Vol. 9, No. 4, p 419-429, October 1972 


The computer program can be used for the preliminary design of 
a propeller having a prescribed pitch dis.ribution. The viscous 
effects of the propeller are taken into account by giving as 
input the blade-section drag coefficient and chord lengths. 
Blade stresses are computed usina simple beam theory by giving 


as additionai input the blade thickness, rake, and skew. 


DESCRIPTION OF INPUT DATA 


Dimensioned sarameters may be input in the appropriate $.1. or 
English units, as specified by the user illustrated in Table 1. 
Effective power, shaft power, speed (v), number of blades (Z), 
diameter (DB), rpm, propeller wake (1-w,), and estimate of the 
hydrodynamic flow angle distribution (8,) are required input 
parameters in order to make nonviscous propeller design 
calculation 1,2, The radial distribution of blade chord lengths 
nondimensionalized on diameter (c/D) and section drag coefficients 
(Cp) must also be specified as input if design calculations are 
to account for the propeller viscous effects. Since the computer 
program presented calculates the propeller principal stresses 
based on beam theory, the radial distribution of maximum thickness 
nondimensionalized on chord lenath (t/c), rake nondimensionalized 
on diameter (Zp/D), and skew angle (80) must also be input parameters. 
The blade section cavitation number (o) is a required parameter in 


order to make blade surface cavitation checks using the Burrill] 


cavitation diagrams of Reference 7 and Brockett's incipient cavitation 


eee ome ee +e 
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diagrams of Reference 8. By specifying the static head (H) as input 
data, the cavitation number (co) is calculated on the computer. 

Hub dimensions can be input or assumed to be a circular cylinder 

of equal length and diameter as described in Appendix A. A brief 
description of how these input parameters can be determined will 


be discussed in this section. 


Effective Power, Speed, and Shaft Power 


Effective power and speed are normally obtained from model 
self-propulsion experiments. Input effective power (Pe) and shaft 
power (P.) are defined as follows: 

Pe = VT(1-t) 
Pe = 2nnQ 
where n = propeller revolutions 
Ps = shaft power, 
Pe = effective power, 
Q = propeller torque, 
T = propeller thrust, 
V = ship speed 
(1-t) = thrust deduction, which may vary with 


diameter and speed. 
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Nondimensional Radial Distance (Xx) 


This is a reference set of eleven nondimensional radia} 
distances x, at which all other distributions, either input 
or calculated by the computer, are defined as existing. In 
general x, = r/R, with the restrictions 

x, = r/R 

X44" R/R=1, 
where re = the distance along the proneller reference line from 


the shaft axis to the ith section, 


a 


propeller hub radius, and 


x 
u 


propeller tip radius. 
Propeller Wake 


The radial distribution of the axial wake (1-w,) which varies 
with propeller diameter is also required input data. The circum- 
ferential mean of the axial velocity distribution (1-w.) is obtained 
from a wake survey without the propeller operating. However, the 
(1-w.) wake distribution must be corrected for the propeller action. 
No completely satisfactory method is presently available to obtain 
this correction, but an approximation of the radial distribution of 


the wake (1-w,) with the propeller uperating is obtained as follows: 


hile 


Wake distribution 
(1-w,) = C(1-w7)(1-w.)J/(1-4,) 1] 
where (1-w.) = radial distribution of the circumferential 
mean wake from wake survey data, 


(1-w,) = effective wake from self propulsion data 


(1-w,) = volume mean wake, [2/(1-x2) Is} (l-w.) xdx, 
R = propeller radius, 

r = sropeller jocal radius 

rh = propeller hub radius, 

x = nondimensional radial distance (r/R), and 
Xp = nondimensional hub radius (r/R). 


The propeller wake distribution may also vary with propeller 


diameter depending on the hull characteristics of the vessel. 


Advance Angle Distribution Input Option 


The advance angle distribution (tang) is normally calculated 
from the equation V(1-w,)/(mnDx) using the propeller wake (1-w,) 
from Equation [1]. For most single screw propeller designs this 
approach gives good performance predictions. In a few instances, 
it may be desired to utilize the computer program presented to 
design and predict the performance of propellers operating inside 
a duct or in the vicinity of another propeller as in the case of © 


tandem or contrarotating propellers where the axial (w,/V) and 


“je 


tangential (w,/Y) velocities induced by these additional sources 
near the propeller plane can be predicted. In this case, an option 
in the computer program allows the input of the estimated tans 
distribution which accounts for the axial (w,/V) and tangential 
(w./V) velocities induced by other sources in the following manner: 
Equation for Estimating tan& If Input Option Is Used 


tang = [(1-w,) #0, /V)I/E(x/,3-(/V)] (2] 


estimated 


where w,/V axial velocity from other sources, 


if 


w,/V = tangential velocity from other sources, 


X 
Ss 


V 


advance ratio based on ship speed, V/(xnD), 


ship speed, and 


D 


propeller diameter. 

It can be seen from Equation [2] that for the case where (w./V) 
and (W/V) values are specified as zero, the advance angle tans 
is calculated in the usual manner when designing single screw 
propellers. If the tang input option is used, the tans values 


are calculated using Equation (°j. For *he normal sipsle screw 


propeller desiqn case, tang is calculated on the computer. 


Hydrodynamic Flow Angle Distribution 


The hydrodynamic flow angle distribution (tans) can be 


specified as input. An option is included so Lerb's optimum 


10 


tan6, distribution® 


can be caiculated by the computer as follows: 
tangy = (tang/ng)[(1-0)/(1-n,)]/? (3] 
where n. = propeller ideal efficiency 
ny = 0.85 is used in the program, and 
tang = advance angle distribution 
Lerbs' optimum tang, distribution usually results in optimum pro- 
peller efficiency. If other factors such as cavitation, strength 


and vibration are considered, the input of an alternate tanB, 


distribution may be desired. 


Static Head 


The static head (H) at the shaft centerline is 
required input. This parameter (H) is defined as HotH.-H where 
H. is the shaft submergence, He is the atmospheric pressure, and 
Hy is the vapor pressure of fluid which is normally smal] compared 
with H, and may be neglected. The static head (H) is used to 
calculate the section cavitation number (co) in Equation [25] and 


the Burrill cavitation number oo 7 of Equatiun [30]. 


Blade Outline and Expanded Area Ratio 


The blade outline (c/D) and expanded area ratio (Ae/Ay) must 


be input for the design. An expanded area ratio (A-/Ay) is 


calculated on the computer according to: 


1] | 


A-/Ay eae c/D dx 
where c/D = nondimensional chord length, and 

Z = number of blades 

If the (Ap/Ap) input value differs from the one indicated 
in equation [4], adjusted (c/D) values used for the design are 
calculated as follows: 

(c/D) 457 (c/0) input Ae/AQ) catc! (Ae/Ag) design 
The final blade outline and expanded area ratio should be chosen 


to give satisfactory propeller strength and cavitation characteristics. 


Blade Thickness to Chord Ratio 


The input of maximum thickness to chord ratio (t/c) 
values allow an estimate of the propeller principal stresses 
(see The Propeller Stress Calculations Using Beam Theory section 


discussed later) based on beam theory® 


to be calculated during tke 
preliminary design stage of the propellers. From a rough estimate 
of the blade outline (c/D) for the final design and an estimate of 
the radial distribution of thickness (t/D) based on fatique 
strength®, the following equation can be used to obtain initial 


(t/c) input values: 


Blade Thickness Ratio: 
t/c = (t/D)/(c/D) 
where t/D = radial distribution of thickness (can be estimated from 


Reference 6). 


12 
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Rake and Skew 


Nondimensional rake (Z,/0) and the skew angles (8,) for a design 
are specified to permit adequate predictions of principal propeller 
stresses using the beam theory method described in Reference 6 
and discussed later in the propeller stress section of this report. 

The rake (Zp) is defined consistent with Reference 1] (see Fiqures 
1 through 2d) as the distance from the propeiler plane to the genera- 
tor line in the direction of the shaft axis. Aft disolacement is 
considered positive rake (see Figures 2a). 

Since the skew angles (8,) significantly affect propeller 
unsteady forces, a computer program based on the unsteady 


propeller lifting surface theery of References 12, 13, 14, and 15 


pee 


can be used to select the skew angles (80) for the propeller 
design. The input skew angles (8,) in degrees are defined as 
the angular displacement of points on the blade reference line 


from the propeller reference line in the projected view. 


11. Cumming, R.A., Dictionary of Ship Hydrodynamics - Propeller 
Section, 14th International Towing Tank onferenc: OTS, 
Report of Presentation Committee, Appendix VII 

12. Tsakonas, S., Breslin, J., and Miller, M., "Correlation and 
Application of an Unsteady Flow Theory for Propeller Forces," 
Transactions of the Society of Naval Architects and Marine 
Engineers, Vol. 75, p 158-193, 1967 

13. Breslin, J.P., "Exciting-Force Nperators for Ship Propellers," 
Journal of Hydronautics, Vol. 5, No. 3, p 85-90, July 1971 

14. Jacobs, W.K., Mercier, J., and Tsakonas, S., "Theory and 
Measurements of the Propeller-Induced Vibratory Field," 
Davidson Laboratory Report SIT-DL-1485, December 1970 

15. Tsakonas, S., Jacobs, W.R., and Ali, M.R., "An Exact Linear 
Lifting Surface Theory for Marine Propeller in a Nonuniform 

i Flow Field," Journal of Ship Research, Vol. 17, No. 4, 


December 1974 
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Section Drag Coefficient 
In order to account for viscous effects when predicting the performance 
of a propeller, the section drag coefficient (Cy) must be specified as input. 
A section drag coefficient (Cp) value of 0.0085 usually gives reasonable 
estimates of model propeller drag for propeller shapes normally used at DTRSRDC 
in the past. For propellers having very thick blades, the following equation, 
available as an input option on the computer, and derived as a function of 
maximum thickness (t/c) values using experimental data from NACA 66 type 


sections, (07 !? 


Section Drag Coefficient: 
4 
Cy = Cry (1+1.25(t/c)+125 (t/c)°] 
where Ce is the friction resistance of the section, e.g., Cro = 0.008 for 


° 


Reynolds number of approximately 10° and Ceo = 0.004 for Reynolds numbers of 


approximetely 10°, 


Options for using alternate nonlinear Ch distributions, or a constant Cy 


distributions are also available. 


Final Piten Ratio (P/D) Input Option 

The final pitch-diameter ratio (P/D) distribution defined as 
nxtan ¢ where » is the pitch angle is normally not known 
during the preliminary design stage of the propeller. Because 
of this hydrodynamic flow angle (8;) from lifting line 


16. abbot, Ira H. and Von Doenhoff, Albert E., “Theory of Wing 


Sections Including a Summary of Airfoil Data," Dover Publications 
Inc., New York, Library of Congress Catalog No. 60-1601, 1949. 


Hoerner, S.F., "Fluid-Dynamic Drag," Published by the author, 
Midland Park, New Jersey, 1965. 
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will give a better estimate of the section drag coefficient (Cp): 


ional 


P wailing, 


i 
‘ 


calculations is substituted for the pitch angle (¢) in the 


bending moment calculations when making stress calculations 


using beam theory and when predicting clearance between blades 


and fillets at the propeller hub. 


known, they can be used. 


If the final P/D values are 


DESIGN AND PROPULSIVE PERFORMANCE PARAMETERS CALCULATED 


In order to simplify propeller design calculation pro- 


cedures, the desian thrust and power of propellers are usually 


considered in nondimensional form, 


The nondimensional design 


thrust loading coefficients (Cre and Cry) and design power 


joading coefficients (Coe and C,) are calculated on the computer 


in the following manner: 


Thrust Loading Coefficients: 
Cog = T/L( 2) IRV" or Coy = T/L (/2)aR?v," [7] 


Power Loading Coefficients: 


Cbs 


where V 
Va 
p = 


ship speed, 


speed of advance of propeller, V(1-wr), and 


density of fluid 


Lerbs lifting line theory is used to determine the 


propeller lift coefficient (C,), circulation (G), hydro- 


dynamic flow angle (8,), axial induced velocity (U,/2V), and 


en ae NLC, Bite en = 


tangential induced velocity (U,/2v). These lifting line calcula- 
tions take into account propeller viscous effects by specifying 
as input in the computer program the propeller section nondimen- 
sional chord length (c/D) and section drag coefficient (Cp). A 
method for obtaining values for c/D and Ch is discussed in the 
Description of Input Data section of the report. Design parameters 
such as thrust loading coefficient (Cro), thrust power coefficient 
(Cri) power loading coefficient (Co) propeller efficiency (np), 
and oropulsive efficiency (np) are calculated for each propeller 
in the following manner: 
Thrust Loading Coefficient: 
252 
=fl « = 2 

Crs ‘xp ( etan8,) (dC p< /dx)dx T/L( 6/2) RV" [9] 

Thrust Power Coefficient: 
fue 
=i 2 ea! = - 

Conp xy w)(1 etan8,) (dC 1. ;/dx)dx (1 wr ){T/[(¢ /2)nR°V°]} [10] 

Power Loading Coefficient: 
=f e 243 

Cose A AI d tang )dCy¢./dx)dx=P./[( 0 /2)nR°V~ J [11] 
Estimated Propeller Efficiency: 

Npe= Copp/Cpe * TVa/(2nQn) [12] 
Estimated Propulsive Efficiency: 


np=(1-t)Cro/Cpo=Pe/Pp [+34 
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where C = section lift coefficient from lifting line theory 
Cro: = nondimensional inviscid local thrust loading coef- 
ficient, arg x/rt,/2V] 
Coss = nondimensional inviscid local power loading coef- 


ficient, (4Z/A)xG[(1-w, )+U7/2V] 
G = nondimensional circulation from lifting line theory 
U,/2v = axial induced velocity from lifting line theory 
Ur/2v = tangential induced velocity from lifting line theory 
The calculated propeller thrust (T) is obtained by substituting 
the Crs calculated in Equation [9] inte Equation [7] and the 
delivered power Ph is obtained by substituting Cos computed 
using Equation [11] into Equation [8]. 
{ Other parameters useful in designing and evaluating the 
| performance of propeliers include the advance coefficient (J), 
ship speed advance coefficient (Jy). thrust coefficient (Ke), 
torque coefficient (Ky) moment due to thrust (Mo) moment 
due to torque (Mop) » moment parallel] to section nose-tail 
line (Mo) moment perpendicular to the nose-tail line (M9) and 
the blade loading distribution (Ly). These parameters are 


calculated as follows: 


Advance Coefficient: 


J=V(T-w7)/(nD)=V,/(nD) [14] 


tae owe o 


Ship Speed Advance Coefficient: 


Jy=V/ (nD) {15} 
Thrust Coefficient: 
KeT/( ¢0°04)=(nCre/8)0y° [16] 
T ayys ty 
Torque Coefficient: 
s 3 
K-=0/( ¢ n2D?)= (Cpg/16)dy [17] 


Moment Due to Thrust: 
Mop (x, )=L e wR VED (22)12 (X-¥5) (1 e tang, )[dC,, «/dx dx [18] 
Moment Due to Torque: 
_ 3,,2 } 
Mop(Xp)=Le mR°Y HAZ) Uy, (x-xq)(tandy+ )[dCy¢ 5 /dx dx [19] 
Moments Parallel to Section Nose-Tail Line: 


Mg (Xg)*Mypcoso+Mo, sing [20} 


Moment Perpendicular to Section nuse-Tail Line: 


Mg (Xp)=ty,sine-My,cose {21] 
Blade Loading Distribution: 
L(x) ~(1/2)ocv,.“C, [22] 


where X5X_ = propeller nondimensional radia} stations, 
r/R and r/R 
$ = propeller pitch angle (input or 81) 
V.= section inflow velocity, VE (1-w, )+Up/2V1#L (x72, )-U,/2V1 
Appendix A presents procedures for calculating propeller weight, 
center of gravity, mass polar uments of inertia and radius of 


gyration and Appendix B shows how the American Bureau of Shipping 


(ABS) coefficients and thickness requirements are calculated. Other 
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calculations made include stress calculations based on beam 
theory, parameters for making blade surface cavitation checks, 
chord lengths for making final pitch and camber calculations, 
and the prediction of spacing between blades and filiets at 


the propeller hub. 
PROPELLER STRESS CALCULATIONS USING BEAM THEORY 


A propeller blade must contain enough material to keep the 
stresses within a blade below a certain predetermined level. 
This level depends on the material properties with regard to 
both steady-state and fatigue strength and to both mean and 
unsteady blade loading. The material selection controls the 
allowable stress level and the blade chord, thickness, rake 
and skew are the main parameters which contro] the blade stress 
for a given blade loading. The principal stresses in the 
propeller blade are computed for each propeller. Both hydrody- 
namic and centrifugal loadings are considered. In this stress 
calculation procedure, the propeller blade is represented as a 
straight cantilever beam of variable cross-section without camber. 
Experimental results show that the neutral axis of an airfoil 
section lies approximately along the mean line so camber is not 
considered in the stress calculations presented. Only the 


maximum principal stresses calculated at the mid-chord of each 


section are printed as output in the computer program. Stresses 


for the final design should be calculated by finite element 
techniques if rake and skew for the propelter differ from 


usual propeller shapes. 
PARAMETERS FOR MAKING BLADE SURFACE CAVITATION CHECKS 


Brockett's theoretically derived incipient cavitation charts 
of Reference 8 can be used to predict the blade surface cavita- 
tion characteristics of each propeller once the lifting line cal- 
culations have been completed. The two-dimensional camber to 
chord ratio (f,/c)s ideal angle of attack in degrees (a,), section 
cavitation number (co), nondimensionalized with the section inflow 
velocity (V.), and the maximum and minimum fluctuating angle of 
attack (osx? Onin) in degrees are parameters that must be deter- 
mined before Brockett's incipient cavitation charts can be used. 
These parameters are calculated as follows: 

Section Maximum Camber to Chord Ratio for NACA a=0.8 Meanline: 

fyy/c=0.0679 C, [23] 
Section Ideal Angle of Attack in Degrees for NACA a=0.8 Meanline: 

aj=1.54 C, [24] 
Section Cavitation Number: 

a=2g(H-xR)/V," [25] 


where g = acceleration due to gravity 


H 


static head at shaft centerline (see section 
on Static Head). 


The maximum and minimum fluctuating angles of attack 


& emer) 


( 


by Lerbs and Rader in Reference 18. These calculations can be 


Cnax* min? in degrees are calculated using the method derived 


made using the following equations: 


Maximum Fluctuating Angles of Attack: 


Omax 04” (48) F(x) [26] 
HNinimum Fluctuating Angles of Attack: 
Onin oy” (+48) F(x) [27] 


where -a8 = Maximum effective angle of attack in decrees 
(from wake survey data), and 
+a8 = Minimum effective angle of attack in degrees 
(from wake survey data). 

The parameter F(x) in Equations [26] and [27] is dependent 
on the hydrodynamic flow angle (87), the advance angle ‘g) and 
the lift coefficient (C)s and is calculated on the computer 
using the following equation: 

F(x)=1/(1+2ntan(B,-8)/C, J (28] 

The Burrill] Cavitation Charts’, can also be used to give an 
approximate check on the cavitation performance of propellers. 
Burrill's thrust loading coefficient (t.) and cavitation number 
(09,7) at 0.7 radius defined as follows are parameters that 
must be known to use these cavitation charts. 


18. Lerbs, H.W. and Rader, H.P., "Uber de Auftriebsgradienten von 


Profilen im Propeller Verband," Schiffstechnik, Vol. 9, No. 48, 
p 178-180, 1962 


2] 


reeT/{(0 /2)ApL(W(T-Wy 4 9) )°4(0.79n0 17) [29] 
Burrill Cavitation Number at 0.7 Radius: 
0, 7=20H/L{V(1-w,_,))°#(0.79nD}7] [30] 


where A. = propeller expanded area, Re dr, 
h 
Ay = propeller disc area, nD°/4, 


A 


p = propeller projected area, (1.067-0.229(P/D; JA... and 


(P/0), 


propeller pitch ratio at 0.7 radius input or 
estimated as 0.?ntans,. 

Keller's method? of predicting the minimum expandea area 
ratio of the propeller is also calculated on the computer. Tne 
minimum expanded area ratio, based on eliminating back bubble 
type cavitation, is computed as follows: 


(Ag/Ag) y=K_(2-6+0.6Z)Kq/ {oq 71I-#(0.78) 73 [33] 


The constant K=0.15 was used in this program. Many possible 
alternate values of K may be desired, for example: K=0.2 is 

used for single screw ships with bronze propellers having rake 

of approximately 10 degrees, K=0.10 is used for twin screw ships 
with copper-aluminum propellers, K=0.15 is used for twin screw 
ships with bronze propellers and for single screw ships with 
copper-nickel-aluminum propellers, and K=0 to 0.05 is used for 
propellers for fast ships such as destroyers and frigates. If 

a different value of K is desired, the expanded area ratio cal- 
culated in equation [31] should be adjusted to account for changes 


in the value of K. 
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CHORD LENGTHS FOR PITCH AND CAMBER CALCULATIONS 


The final pitch and camber for each propeller can be calcu- 
lated using computer programs presented in References 3, 19, 20, 
21, 22, and 23. The programs require as input the section chord 
length, (c/R), - and (c/R) rp nondimensionalized on propeller radius, 
in terms of the skew angle (8), hydrodynamic flow angle (81) and 
blade outline (c/D). The parameters (c/R) and (c/R) rz, measured 
from the leading edge and trailing edge to its reference line, 
respectively, are calculated in the following manner on the 


computer: 


19. Kerwin, J.E., "The Solution of Propeller Lifting-Surface 
Problems by Vortex Lattice Methods," Department of the 
Naval Architecture and Marine Engineerina, Massachusetts 
Institute of Technology, Cambridge, Massachusetts, 1961 

20. Kerwin, J.E., and Leopold, R., “Propeller Incidence Cor- 
ection Due to Blade Thickness," Journal of Ship Research, 
Vol.7, No. 2, 1963 

21. Cheng, H.M., “Hydordynamic Aspects of Propeller Design Based 
on Lifting-Surface Theory: Part I - Uniform Chordwise Load 
Distribution," David Taylor Model Basin Report 1802, 1964 

22. Cheng, H.M., “Hydrodynamic Aspects of Propeller Design Based 
on Lifting Surface Theory: Part II - Arbitrary Chordwise 
Load Distribution, " David Taylor Model Basin Report 1803, 
1965 

23. Kerwin, J.E., “Computer Techniques for Propeller Blade Section 
Design," Transactions of the Second Lips Propeller Sy-posium, 
Drunen, Holland, p 1-31, May 1973 
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Chord Lengths Measured From Blade Leading Edge: 
(c/R), -=x8/ (57.296 cos¢)-c/D [32] 
Chord Lengths Measured From Blade Trailing Edge: 
(c/R)7¢=X8c/ (57.296 cos$)+c/0 [33] 
(98, is used when pitch is not input) 


SPACING BETWEEN BLADES AND FILLETS 
Propeller designs should have enough clearance between 
blades at the hub so fillets are properly applied. Hill, 
Reference 24, derived the following equation, which is used 
in the program to estimate spacing between blades at the hub 
without fillets: 
G, = (2ur,)/Z - (t,/sing) [34] 
where ¢ = input pitch or By 
Based on a number of full-scale propellers built with standard 
fillets, Hill's blade clearance equation was modified in the 
following manner to estimate spacing between fillets at the hub 
during the preliminary stage of the design: 
G. =(énr,/Z) - (1.9t,/sing) [35] 
A layout of blade sections 1s recommended as a final fillet 


clearance check. 


PROPELLER DESIGN THRUST AND POWER OPTIONS 


The thrust option can be used to make lifting line calcula- 


tions for propellers required to produce a given thrust at 


meine ae ren RY PS 
ww 


24. Hill, J.G., "The Design of Propellers," Transactions of the Society oe 
of Naval Architects and Marine Engineers, Vol. 57, p 143-170, 1949. 
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specified values of speed and rpm (this is accomplished by adjust- 
ing the input tan By distribution), or the power option can be used 
if the propeller is required to absorb a specified power at a given 
rpm (in which case the speed is determined). 

From each calculated power, a new value of speed (assumed to 
vary as the cube root of the ratio of the design and calculated 
power) is obtained and its correspondinn effective power is cbtain- 
ed from the effective power input curve. Design calculations again 
produce a new calculated power, and the process continues until the 
closeness criteria of design calculated power is satisfied (two 
iterations are normally sufficient). Smaller increments of input 
speeds in general cause faster convergence. 

Once the basic shape of the distribution is defined (see the 
Hydrodynamic Flow Angle Distribution section) the final tans, 
distribution Ky ‘anB, is determined using the thrust or power 
options, making lifting line calculations of three nondimensional 
thrust loading coefficients (Clo)y> (Cro)o» and (C1); that cor- 
respond to three hydrodynamic pitch distributions, Ky tang,, Ko 


tangy, and Ky tans, where K = 0.975, Ko = 1.0, and K,= 1.025. 


3 
Once these calculations are obtained, the following system of 


equations are set up: 


7 2 
z 2 
_ 2 
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from which values of A, B, and C are obtained. Then, values of A, B, 
and C are substituted in the following equation to obtain the value of 


K 


4° 
= 2 
Cr. = A+BK,+CK, 

vihere C__ = design thrust loading coefficient (Equation (7)). 


Ts 


COMPUTER PROGRAM AND SAMPLE DESIGN 
CALCULATIONS 

Liftina line theory of References 1 and 2 has been used to 
derive a computer program for the areliminary design of propellers 
using high speed CDC computers at DTNSRDC. A core size of approxi- 
mateiy 60,000 octal is required for the program. The average run- 
ning time for a design based on thrust is 25 seconds and 100 seconds 
are required using the power option. As mentioned earlier, Appendix 
A presents procedures used to calculate propeller weight, center of 
gravity, mass polar moment of inertia, and radius of gyration. 
Appendix B shows how the American 8ureau of Shipping (ABS) coefficients 
and thickness requirements are calculated. A detailed description of 
the input and output formulas for the computer program is present in 
Appendix C exercising most of the available options. Design calculations 
are presented in Appendix D for a sample design using the thrust 
option and results using the power option are shown in Appendix E. 


A Fortran listing of the computer program is presented in Appendix F. 


26 


V, (9) 


{wr = Vo (r,0)] 


Figure | — Typical Velocity Diagram for a Propeller Blade Section at Radius r 


(The diagram is drawn wit. all quantities positive and the velocity vectors represent the velocity of 
the propelicr blade ection relative to the Muid) 
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Figure 2b - Diagram Showing Recommended Reference Lines (Looking to Part) 
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PROPELLER REFERENCE LINE 
AND GENERATOR LINE 


BLADE REFERENCE LINE 
(LOCUS OF BLADE SECTION 
REFERENCE POINTS) 


PROJECTED BLADE 
OUTLINE 


TRAILING EDGE 


cba 
a 


LEADING EDGE 


PROPELLER HUB REFERENCE POINT OF ROOT SECTION 


STARBOARD. 
DOWN 


= 


SHAFT AXIS 


NOTE: THE SKEW ANGLE, @,, SHOWN AT RADIUS r IS LESS THAN ZERO. 


Figure 2c = Diagram Showing Recommended Reference Lines (Looking Forward) 
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FIGURE 2d - INPUT FOR DETAILED HUB OPTION 
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TABLE 1 
List of dimensioned input and output parameters used by 


computer program based on English and SI units 


(1) 


Parameter English Units SI Units KSI 
Shaft power (Ps) hp KW 0.7457 
Effective power (P,) hp KW 0.7457 
°, Lbm/£t? kg/m" 16.01846 

(2) s 
Vv knots n/s 0.514444 

(2) 

Vv ft/sec knots 0.5924 
D ft m 0.3048 
H ft .) 0.3048 
fe) lbf Sec" /Eee kg/m? 515.3788 

‘ (3) 
n rev/min rev/min 1.0 
T, weight lbf N 4.44822 
ve ft/sec m/s 0.3048 
LI lbf/f£t N/m 14.5939 
Mop in lbf Nn 0.112585 
Mop in lbf Nn 0.112985 
M in lbf Nn 0.112985 
xO 
M in lbf Nn 0.122985 
yo 
Max Stress bE/in® Pa 6894.757 
SKEW deg deg‘) 1.0 
RAKE deg deg >? 1.0 
Mass polar moment lbm in® kgm 0.00029264 


of INERTIA 


(1) Multiply English Units by KSI to get SI Units. 

(2) Computer program uses knots in both English Units option 
and SI Units option. 

{2) These are not SI Units but are permitted to be used in the 
SI syet*™ according to International Standards Organization (150) 
Standard No. ita. 
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APPENDIX A 
PROPELLER WEIGHT, MASS POLAR MOMENTS OF INERTIA, AND CENTER 
OF GRAVITY 


The approximate propeller weight (HW) and location of center 
of gravity (CG) from the propeller center line is also calculated 
for each design. In order to make these calculations, the density 
of the material ( ) must be specified as input. 

Hub dimensions can be input, or a solid circular cylinder of 
equal length and diameter, with propeller center line at the mid 
length of the hub can be assumed. If the cylindrical hub is assumed: 

No = Wath, 
where W, = weight of blades, ply a(x)dx 


att 3 
W, = weight of hub, (nopD,/4), 
Di, = hub diameter 


a(x) = area of section at radius x, 2c(x)t(x) s ! t(x,x, )dx, 
c(x) = chord length of section at radius x, 
t(x) = maximum thickness of section 


t(x4x,) = Chordwise distribution of section thickness (NACA 66 
modified thickness form is used in program) 


Xo nondimensional coordinate alona the section chord, (r/c), 
and 
= density of material considered 
33 : 


w 
wm 
~ 


3 pee ty SE SS oy hi 
oe 


if calculations are desired for a noncylindrical hub, the follow- 
ing dimensions are input: 


Dey = forward hub diameter, 


© 
iT] 


A aft hub diameter, 


- 
" 


length of the hub, 


Dep = forward bore diameter, 


rc 
i] 


R distance from aft end of hub to propeller reference line 
The hub is then treated as consisting of two frustrums of cones, 


joined at the propeller reference line. Again 


Wo = Wp +H 
but now 
Wy = Waptle eM 


2 
where Wye = weight of solid aft frustum, (np 2/3) (Dy,,0,,+D,,+0,,)/4 


Nee = weight of solid forward frustum, 
2 2 


H H FH 
; 2 2 
b 7 weight of bore, (n¢ p/3) (Dip tOnpDeptOee )/4 


The propeller center of gravity (CG) with respect to the blade 


= 
it 


center line, where plus values represent the distance ahead of the 
center line and negative values aft of the center line, is calculated 
in the following manner (see Figures 2b and 2c): 
Center of Gravity: 

CG = Mf M 


where M = moment of the propeller, cell, a(x) B(x) dx 
h 


grain 


1 


where M, = moment of the propeller, %p qs a(x) B(x) dx 
h 


B(x) = distance of center of gravity from propeller reference 


line, x +y +H 
x = the distance of the center of gravity of the section along 


the chord line from the position of maximum thickness, 


y = the distance between the center of gravity of the section 
and the chord line of the section, measured perpendicular 
to the chord line, 


r. = the radial distance from the shaft axis to the section, 


R= rake angle in radians of the section, 
OF skew angle in radians of the section, and 


¢ = pitch angle in radians of the section ( ¢=8, is used 
when pitch is not input). 


x<| 
" 


Then the component of the center of gravity adjusted for 
pitch, xsin¢ +ycos ¢ ; 
y = the component of the center of aravity adjusted for skew 


(=- Gritan¢) and rake (=- @ r./R); and 


H = the component of the center of gravity along the shaft 
axis of the specified hub, aft of the propeller reference 
line taken as positive. 


35 


7; 
.¢ 


aa erences ee, 


If the standard hub is selected, 
= a 2 
hon = mass polar moment of hub = WO, /8 
otherwise, 


H Agat F 


po = ™ © p(Apis*FpynBpyy) /10 
where Aow = aft mass polar moment = 


43 2 3. 4 
Lay (04, +0," 09).+D, Dani? Pay +0 yy, )/16 


Fone = forward mass polar moment = 


440, 3p_,.40,2D-.,-+#D,,D 


Dy PO Oey Oy, Dey FO Dey, 


3 
a) ( H 


4 
(Hj -L +Dey')/16 


Boy = mass polar moment of bore = 


4. 35 +D 2 2 D 


3 
Bret ap ea + ap? Fe 


4 
H D +0cp )/16 


Lvs 


Zon = mass polar moment of blades = 


See 
“4 ok i? € (x) dx 
Then Pow = propeiler mass polar moment= 


Zon Hon 


and Ks = radius of gyration of blades = 


VZou/Wp 


at 
PM 


Ky = Radius of gyration of hub = 
V Hoy] ay 
Kp = radius of gyration of propeller = 


¥ Pon! Wp 


why 
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APPENDIX B 
AMERICAN BUREAU OF SHIPPING (ABS) COEFFICIENT AND THICKNESS 
REQUIREMENTS 
The minimum blade thickness at the one-quarter radius is 
defined in Reference 25 by: 
ty=K, YAP./(C, CZn) + C.BZ,/(4C,,C) 
and T3=K, YAP /(C C20) + C.BZ7/(4C,C) 


where Pe = shaft power at the maximum continuous rating, 


n = propeller revolutions per minute rpm at the maximum 
continuous rating, 


Z = number of blades, 
Zp = rake, positive aft, at the blade tip. 


{ Zp i= rake and rake plus skew induced rake, at the blade 
: tip, positive aft taken as positive, 


5 = a/(Woty o6) = section area coefficient at the one- 
quarter radius 


2 
0.25) = section modulus coefficient at the one- 


quarter radius 


C= 1 /(UMet 


A =1+(K,/tan ¢ 0.7)*kgtan 4 0.25 


B =(Ky(Ag/Ay)/(Z)  (n/190)* (0/20)° 
t =(1+Ke tan q 0.25) WoF-B) 


where Ky through Ke are constants which depend on the system of 


units employed 


in 25. American Bureau of Shipping, "Rules for Building and Classing 
( Steel Vessels," p 621, 1972. 
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a, = area of expanded cylindrical section at 0.25R = 


2c t(0.25,x) dx, 


A 
0.259.25 “0 


We = expanded width of cylindrical section at the one-quarter 
radius, 


0.25 = maximum thickness at the one-quarter radius, 


— 
i" 


0 moment of inertia of the exnanded cylindrical section at 
one-quarter radius about a straight line through the center 
of gravity parallel to the chord line, 


Cc 
" 


F maximum normal distance from the moment of inertia axis to 
points on the face of the section, 


Co 25 = section chord length at 0.25R, 


Gaon $9.77 the pitch angles at the one-quarter and seven- 
tenths radius (or Br if ¢not input), 


t(0.25,x) = chordwise distribution of section thicknesses at the 
0.25 radius 

A-/Aq = expanded blade area divided by the disc area, 

D 


propeller diameter, 


f = material constants in length/weight as a function of w 


smal og 


as illustrated in the table below. 


Materials 

Manganese bronze 
Nickel-manganese bronze 
Nickel-aluminum bronze 
fn-Ni-Al bronze 


Cast iron 
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a nondimensional parameter associated with the material type 


Units 
W 
0.30 
0.29 
0.27 
9.27 
0.26 


If Cy # 0.1, the thicknesses t, and t, are recomputed with 


 Recpaten tout ead el sibardih? cate oe ORI 
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APPENDIX C 


DESCRIPTION OF INPUT AND OUTPUT 


Data are input to the program in unformatted groups. Within 
each group the data must be in the order specified, and each 
value must be separated from al] others by a blank or comma. 
Variable names beginning with I,J,K,L,M, or N require an integer 
input, but other variables which happen to have integer values 
do not require a decimal point. Zeroes must actually be keypunched. 
Beyond this there are no restrictions on how many or few data 
cards are needed for a group, except that each group must start 


with a new card. The input required is as follows: 
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Group Parameters Description of Input 
1 10D Number of data sets(each of the groups 
2 through end form a data set) 
2 UI ,U0, Title Use SI in columns 2,3 to pecify input 
in SI units. Use SI in columns 4, 5 
for output in SI units. The remaining 
66 columns are available for identifi- 
cation of input data. 
3 SHPLFL/T] Nesian shaft power or 9 if the thrust * 
option is desired 
TANBI Use 1 unless use of an exact (not ad- 
justed by program) tang, distribution 
is required, then use 0 (see hydrodynamic 
flow angle distribution 
TANG Use 1 to input tan¢ distribution, or 0 
for machine calculations according to (2) 
XPS Use -] to input a skew distribution, or 
input the skew at the tip for a linear 
skew distribution (both in degrees) 
RAKE Use <RAKE< 0.01 to input a rake/) * 
distribution or input the rake/D at 
the tip for a linear distribution 
POO Use PDO 0 to input a P/D distribution, 
otherwise n xtang, is used as an approxima- 
tion 
* As used within the Fortran Program 


4) 


4 
4 
Group 
3 
tran eee, ae 
4 


Parameter 


CD 


TYPE 


HUB 


DIAMEL] 
EWAME 
ETHRUS 


Description of Input 
Use CD>10 to input the radial distribution of 
drag coefficients (CD); 0<CD<10 to input a con- 
stant drag Cp)=CD at all radial stations; CD=0 
causes the computer program to calculate the 
radial distribution of drag coefficients using 
the equation Cy°0.008[141..25(t/c)#125(t/c) "Js 
-10<CD<0 causes the computer program to use a 
constant frictional resistance Ces ABS (Ce) in 
Equation 6 Co-CeoL I+1.25(t/c)#125(t/c) "J; 
CD<-10 to input the radial distribution of 
frictional resistance (CFO) values to be used 
in Equation 6. 
For ABS minimum blade thickness calculations 
coefficients use the number corresponding to 
the material type in the blade thickness 
section. Use TYPE=0 to suppress these calcu- 
lations 
Use HUBAO to input actual dimensions of HUS, 
or HUB=0 to use solid cylindrical hub in 
relevant calculations 
Propeller diameter 
Effective wake (1-7) 
Thrust deduction (1-t) 
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Group 


Parameters 


HEAD[L] 


DEN(M/L3] 
RHO[M/L3] 
IVV 
VEI(L/T] 
EHP[FL/T, 


12Z 
JBL 
IEA 
EXX 


IRPM 
XMM 


ae Nar Ante fe ing Me 


Description of Input 
See section on Static Head, at the shaft 
centerline 
Density of propeller materials 
Density to be used for water 
Number of velocities input 
TYV velocities 
IVV design effective power values, corres- 
pondina to the input velocities 
Number of different blade numbers input 
1ZZ values of blade numbers 
Number of expanded area ratios input 
IEA values of Ac/A,. A value of 0 input as 
the first A-/A, will result in calculations 
being done exactly at the design AL/A,- If 
this default is used, calculations will not 
be done at other values of A-/A, in the 
list that happen to be within 0.005 of the 
design A-/A.. 
Number of different rpms input 


IRPM values of revolutions per minute 
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Parameters 


Description of Input 


(The remaining five radial distributions have eleven values each) 


X3 


x4 


X5 


X6 
AZZ(25) 


Values of r/R at which other radial 
distributions are defined 

Propeller wake (1-w,) 

Hydrodynamic fiow angles (tane,). If zeroes 
are input, Lerbs optimum is calculated as a 
default distribution 

Section chord lengths (c/D) 


Thickness to chord ratios 


(Groups six through ten which are specified in group three are 


defined at the eleven radial stations) 


AZZ(24) 


FWDDIAM[L J 


AFTDIAM[L] 
HUBLEN(L] 
FDBORE(L] 


ADBORE[L] 
DISREFL[L] 


Blade skew angles, degrees (if XPS<0) 
Section drag coefficients (if CF>10) 
or frictional resistance of section 
Ceo (if CD<+10 

Tangents of advance angles (if TANB>O) 
Rake/D, aft positive (if 0<RAKE<0.01) 
Pitch, P/D (if PD>0) 

Dimensions, as indicated in Figure ed 


(if HUB#0) 


4 
44 3 


A Description of the Output Generated by the Program Follows: 


Program Output Description 

G Nondimensional circulation 

UT/2V Tangential velocity induced at lifting line 
UA/2V Axial velocity induced at lifting line 

DCTSI Local nonviscous thrust coefficient (Defined in 


Equation (9)) 
DCPSI Local nonviscous power coefficient (Defined in 


Equation (11)) 


VR [L/T] Section inflow velocity, equation (22) 
CAVV Section cavitation number, Equation (25) 
CPTI tlonviscous thrust power coefficient (Equation (10) 
when ¢€30). rhs 
CPSI Nonviscous power coefficient (Equation (11) 
when ¢ =0) 
ETAI Estimated nonviscous propeller efficiency 


(Equation (12) when « =9) 


CTSI Nonviscous thrust coefficient (Equation (9) when 
e =0) 
CPT Thrust power loading coefficient (Equation (10)) 
CPS Power loading coefficient (Equation (11)). 
TETS Projected skew angle in degrees 
RAKG Rake/Diameter 
PE Effective power 
ry 
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Program Qutput 
PS 

ETA 

CTS 

CL 

ALI 


FM/C 


CD/CL 
F(X) 


LI[F/L] 
(C/R)DLE 


(C/R)DTE 


T/RD 


PC 


PS[FL/T] 


DESIGN [F] 


CALCULATED [F] 


Description 


Shaft power 


Estimated propeller efficiency (Equation (12) 
Thrust loading coefficient (Equation (9)) 

Section lift coefficient 

Sectien two-dimensional ideal angle of attack 

in degrees for NACA a=9.8 meanline, Equation 

(24) 

Section two-dimensional maximum camber ratio 

for NACA a=0.8 meanline, Equation (23) 

Section drag-lift ratio (e=C,/C, ) 

Parameter for calculations section fluctuating 
angles of attack, Equation (28) 

Propeller blade loading distribution Equation (22) 
Chord lengths for lifting surface pitch and 
camber calculations, Equation (32) 

Chord lengths for lifting surface pitch and 
camber calculations, Equation (33) 

Ratio of section thickness to radius 

Estimated propulsive efficiency, Equation (13) 
Calculated shaft power delivered at the propeller, 
Equation (13) 


Design thrust 


Calculated thrust, Equation (7) and (9) 


2 eat ares 


Program Output Description 


AEPA 17] Area of section 

XBAR [L] Longitudinal position about x axis parallel 
to nose-tai] line from centroid 

YBAR [L] Vertical distance about y axis perpendicular 
to nose-tail line from centroid 

1X0 14] Moment of inertia about y axis perpendicular 


to nose-tail line 


MXO [FL] Bending moment about the x axis, Equation [20] 

MYO [FL] Bending moment about the y axis, Equation [21] 

MTB [FL] Bending moment due to thrust, Equation [18] 

MQB [FL] Bending moment due to torque, Equation [19] 

MAX STRESS (F/L2] Maximum stress 


WEIGHT OF BLADES [F] Equation [32], W,=0 
WEIGHT OF PROP (BLADES 

+DESIGNATED HUB) [F] 

(CENTER OF GRAVITY OF Appendix A 

PROP/D) 

(CENTER CF GRAVITY OF Appendix A 


BLADES/D) 
KELLERS MINIMUM EAR Equation [31] 
SPEED COEFF(JS) Equation [15] 
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Program Output Description 


Advance Coeff (JA) Equation [14] 


THRUST COEFF (KT) Equation [16] 


TORQUE EOEFF (KQ) Equation [17] 


PROPULSIVE EFFICIENCY Equation [13] 
(PC) 


or THRUST COEFF Equation [29] 
TC) 


BURRILL CAVITATION Equation [30] 
COEFF 


CLEARANCE AT HUB BETWEEN Equation [34] 
BLADES/D 


CLEARANCE AT HUB 
BETWEEN FILLETS/D 


Equation [35] 
MASS POLAR MOMENT OF 2 
INERTIA OF CLADES [FL“] Appendix A 


TOTAL MASS POLAR 9 
MOMENT OF INERTIA [FL™] : 


RADIUS OF GYRATION OF : 
BLADES/D 


RADIUS OF GYRATION OF ‘ 
HUG/D 


TOTAL RADIUS OF GYRATION/D " 


ABS MINIMUM THICKNESS Appendix B 
CONVENTIONAL RAKE 


CONVENTIONAL + SKEW 4 
INDUCED RAKE/D 


VALUES USED IN 
DETERMINING THICKNESS " 


. 


ee a ee ears 


Program Qutput Description 


SECTION AREA COEFFICIENT Appendix B 


SECTION MODULUS ‘. 
COEFFICIENT 


AREA OF EXPANDED 9 
CYLINDRICAL SECTION [L“] : 


49 


APPENDIX D 


SAMPLE DESIGN USING THE THRUST OPTION 


OPTIONS EXERCISED IN APPENDIX D 


Power or thrust 

Calculations with input tang, 

Tang, (Lerbs or input) 

Tang (Calculated from (1-w,) or input) 
Cy (constant, variable, or calculated) 
Multiple RPM 

Multiple Z 

Multiple AE/A0 

Check AE/AO of input C/D and modify 
Skew (Linear/Non] inear) 

Rake (Linear/Nonl inear) 

P/D (input or approximated by mxtang, ) 
ABS coefficients 

Hub geometry 

Input Units 

Output Units 
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OPTIONS EXERCISED IN APPENDIX E 


Power or thrust 

Calculations with input tang, 

TanBy (Lerbs or Input) 

Tang (calculated from (1-w,) or input) 
Cy (constant, variable, or calculated) 
Multiple RPM 

Multiple Z 

Multiple AE/A0 

Check AE/AO of input C/D and modify 
SKEW (Linear/nonlinear) 

RAKE (Linear/non? inear ) 

P/D (input or approximated by mxtanB,) 
ABS Coefficients 


Hub Geometry 


Input Units 
Output Units 


Power 


N/A 
Lerbs 
Calc 
Calc 
Yes 
Yes 
Yes 
Yes 
NL 


NL 


Approx 
Calc 


Yes 


Approx 
Calc 
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FORTRAN LISTING OF COMPUTER PROGRAM 
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PROGRAM GMAINCINPUT=256, OUTPUT=512, TAPFS=INOUT, TAPE G=OUTPUT) 
OIMENSTION CHIR0(11) sTHICKNS(119 ,CAMBER (11) ePITCH(11).SKEWR (11) 
1 oxC1L1L9,-8T C11) 

DIMENSTON 6(364.38@) ,DENS(6? .82(161).83 (181) 
1 ,811(181),81% (181) 
2 eH1S (LAL oS AZ(CLLeiL) eBHUL1,11) 0812517) CH (12,12) 

DIMENSION X3CL19 XGC11) XSCLL De MOCL1) sVEL (9), EHP (9) gSLA(9) pcXX(9) 
OIMENSION AZZ (14438) -ASHP (9) »XMM(9) ,-CAVI9) ,CAE (9) -FX(11),883 (11) 
OIMENSION VEL 109) -EHPL (9) .AX 011) pRAK(11) ,OXTST C119 PXT9011)9,SK07) 

OIMENSION VSUSRSQC1L1) ,VSUI9 (11) 

OIMENSION ARZA (7) -XBAR (7), YSA9(7) ,AVYEXOI7) -AYEVOC7TIVEMKXO(7),=NYO(7 
1).°MTB(75 ,£EN9B(7) »STRMAX (7) 

DIMENSION HUL3) eEMMXCL3) eVTXCL 3) TXOLL) CKCLID SFEMKCL1) JOEL E1L1). 
1 BROMI(CLL) sPMOX(11) 9 F(23DSCL39,P C13) 

OTMENSION TABS(E) 

OIMENCTON O4V 09) OME (9) » 04X09) 

COMMON /UNITS/ SI,UI,U0 

COMMON/CWEIGIT /X% «CHORD» THICKNS »e CAYRER gPT TCH, SKEW? DIAM, ZZ DE Ny CAKE 
1ePIyPP7,PPBy2PG,PP1LEWAKE pVS 9 RPS eSIGMA EAR, ATP PIO 
2 oF WONT AM, A> TOIAM,HUBLEN,FOBORE « AODBORE NISREFL 

COMMGN 49, 82,83, Ri 94144915 ,A7,5H, 
ic.cc. 13 ¢JBsJCeJO,IND ES eCLICis) 

COMMON CCONE(12£) eCCTHO(11) ,CCTHR(11) ,CCFOQ(L1) 

COMMON PP1,P°2,PP3,PP4,PPS5,PPH,PPLG 

COMMON SN(73),.C0(73) 

DIMENSION XR0 20), TANBETI (70) .G(200 e XSL (707 XST (21) WAKE (20) 
1 eUVELA (25) eTHICKP (20) 

OIMENSION TLES (10-2) 

PEAL NeNIe INPUTSLERBSeLINR, NLINR 

DAT& ¢APU,ELRE / 8H DESIGN » 8H INPUT / 

DATL W /10 eg de Becton Be ere e des ve Bog rene rlel 

DATA FAMX 4.0 40 27120 oH B2 905993 5 GES Eo eo F515 Lng 0p S796 g 6 8X92 727, 
1 e 35Abee L713, 6 c/ 
NATE TX Fiery 206669 029075 0 ne HO37 0 0 O95, 6H GH2 oe oe HAST, UOTE Le SIL, 
01877461143, 2.33337 
DATA DENS / 169521. 85651049 Boelerds / 
VATA TABS SLCHABS MINIMU ,134M THICKNESS ,L IHS TN LNOHE ' 


tf 


2 12HS (USING P 10H 9iCH / 
DATLTHS, UBS, VBS»XBS / LCH/D=PTXTANB s1NHT)- ’ 
2 1GH/0 INPUT) - ,16H / 


DATA SHPP,SHPT1,SHPT2,DEN1eDEN2 MENS 0 259 TAI 4 POMS of OMMA, ANY J 
1 /SH SHP=,7H THRUST,7H OPTIONe1 OH, DENSTTY »19HOF POOPI(LE, 
2 BHM/FT3)= 43H 2=_45H EAR=_5H PPM=y2Hy wiTH ADJUSTEN / 


DATA TLES / i0H Xx 0iGH 1-wWyY 9i3H C70 ° 
2 1CH T/¢ eif# TANBI 2iQH TANA ’ 
3 10H TETS(OES) s10H RAKG/N ,10H P/N ’ 
& 1CH CO e AG*1H = / 

DATA INPUT sLERBSsLINR gNLINQSPXTNST,CONST,CALE 7104 YNPUT ’ 
2 1CH LERBS oifK LINEAS si GHNONLIN' AS ) 
3 1CHPIYTANBI es LOHCONSTANT sLOHCALSULATED / 


DATA SIT /2HST/ 

DATA YVE / 7THFT/SEC) / 

DATA EXA,URPM / SHAE/AD »SHIPEV/HMIN) / 
OATA VET eUEF,USL /GH(IN) s4HCFT) » SHIM) = 7/ 
QOATA VL,PE 7 2HV( 4 SHPE( / 

DATA ULEsULS / BHILBF/FT) 4 SHCN/M) / 
DATA USV,UEV / GHM/SEC) ,»6HKNOTS) / 


¢ DSTA USP,UEP / 3HKW) » 3HHP) / 
& CATA UED.USO / 1GH(SLUG/FTS) ,10H(KG/M3) / 
OATA ULES,ULSS »ULE2 pULS2 UNM ,UIL , UPA, UPI N2 J5SHCING) SHUR) 
> SHUINZE»SH(M2) .8H(N-M} sddCIN-LBF) »9H(PA) »9H(LBF/IN2) 7# 
JATA UFE,UFS / SHILBF) ,SHUN) / 
DATA 1192,U03 /1QHOF PROP(KG,BH/M3)= 8 / 
NATTA PHL,U07 /LOH(LBM-IN2)= ,10H(KG—M2)= / 
QLETA G2L,.GR2 / GHRADIUS OF ,.9H GYRATION / 
MATA EFT,U0S / GHFT.=,4HH.= / 
OSTA ELN+J05 /SHINCH) 5 5HM) / 
JATA EINS,UQO /7HINCHES= ,7HMETERS= / 
O8TS £av,SIV /oH RP4= ,6HRAD/S= / 
JATA PINS,JOS / 7HINCHES) »7HMETERS) / 
NATA PINS ,»U09 / 6GHINSHES » GHMETERS / 
DATA EFP,U10 /GHFT)= ,»SHND= / 
7 L122./.025% 
EL221./ELIZELI 
T_F21./.3048 
NoT=1./16.91846 
RHOSI=1./515. 3748 
VST EL. 5 146448 
PWP=1390./745.6939 
SIM=1./.1129848 
UWT=1./6o44822 
JC=11 
NKTS=1.6878 
PIT=3,146159265 
{ ¥12=2.°PI 
32 $1 121,73 
A225,.°(I-1)92.°%°1/369, 
aNOT) =SIN CAC) 
Zi C3CIb=SOS CAC) 
REAUI5S,°9 190 
90 17.72 NNIA=1,190 
S?4=6 ./P12 
READ(5,10908) UI,U0 
Rr AN'S.*) SHP,TANBI» TAND yXPSsRAXKE »2005C0, TYPE sHUB 
O:A0(5,°! OLAM,EWAKE,E THRUS, HEAD, DEN, RHO 
PLATS, 9) TVW, CVELCT) »TaL,IVV) pCEHPO(I) Tai, IVV),IZZ, 
L CJFLGT) »THLe1ZZ) TEA, CEXKCT),T=1,TEA) ,IROM, (XMM(T) ,IT=L,1 RPM) 
ReAD(S oO) CXICLD LHL IC) oCXGCI De LALe JU) CXSTL) pLHL, IC) (XGUIs 
1 IT=1,JC¥,CAZZ(1.25),1T=1,5C) 
IFQXPS.LT Ie) READISs*) CAZZE1526) pL19JC) 
IF (AQS(5D).GE.10.) READIS,*) (9(I,7),I=1,JC) 
17 CTANBOGT.0.) READ(5s°) €3(T.8) Tat, UC) 
IF (DAKE .GT 2.0.6 ANDeRAKE LE 01) READ(5,*) (RAK(ID,I=49IC) 
IF (PMU.GT.2.) READ(S,*) (PUT), T=1,J50) 
I= CHUG.NE 0.) READ(S,*) FWODIAM, AF TOIAM,HUSLENs»FOBORE,ADBORE, 
1 NISPEFL 


C ND OATA ZEAD ANYWHERe BEYOND THIS STATEMENT 


Qrenram 

IFC(UI.NE.SE) PFSDIAN* 12. 
q IF (CRAKE.GT.0. ANDRAKE.LE..01)? GO TO 140 
; RAKE =2AKE *RF 
ci GO TO 142 


110 


140 
141 
142 


241 3 


143 


145 


00 141 T=1eJSC 
RAKCLI=RAK(T) FRE 

UDEN=UED 

IOEN=O0EN 

IFCIOENLT.2? IDEN=4 
IFCIDENSLT.7) DEN=OENSCIOEN) 
UL =UEF 

VUSUEV 

PU=UEP 

IFCUI.NEeSI) GO TO 2413 
IFCIDENeLT.7) OEN=OEN/OST 
UQEN=USD 

UL=USL 

VUSUSV 

PU=USP 

WRITE (6512046) 100 


IE (SHP.£0.0.) WEITE(6,10020) SHPT1,SHPT2 gM=ZN1,0°N2,0EN3,0EN 


LF (SHP.NE*0.) WRITE (6910021) PUs SHP gNEN1,D0EN2eDEN3,DEN 


WRITE (6419927) VLoVUs (VEL(ID T=1 eI VV) 
WRITE (6410045) PEyPU, (EHP(IT) s I=1,I VV) 


WRITE (6519028) UL» OTAMyE WAKE » > THOUS, UL pHTAD, UNEN, RHO 


VUSUEV 

PUSUEP 

UNEN=UE9 

UL=JEF 

ESI=EIN 

FSI=£FP 

ANV=E AV 

SNIP=FINS 

SNTR=FINS 

IF(UO.NE*SI) GO TO 143 
VUSUSV 

PUZUSP 

UNE N=USD 

UL=JSL 

SNIF=UC9 

SNIP=U08 

ANV=SIV 

DEN2=U02 

DEN3=U33 

FSI=uU1C 

ESIT=UC5 

IF(UI.NEeSI) GO TO 149 
OPTION FOR SYSTEM INTERNATIONAL INPUT YNITS 
SHP=SHP*PHR 
OIAM=CIAM*ELF 

HE AD=HEAD*ELF 
RHO=RHOFRHOS I 

00 145 [=isl VV 

VEL CID=VEL(LDFVST 
EHP(LY=EHP(T) *#PHR 
IF(HUB.EO.0.) GO TO 146 
FWOODLAM=FWODITAM*ELF 
AFTDIAM=AFTD TAM*ELF 
HUBLEN=HUBLEN*ELF 
FOBORE=FOBORE* ELF 


wm, 


reaper 


{ 146 


147 
148 
149 


43159 


10645 


26 
7 
10.63 


2 


ADBORE = ADBORE *ELF 
OISPEFLZDISREFL*ELF 


LF (RAKE sGTeQeeANDeRAKE ol Fo201) GO TO 147 


RAKE=RAKE*ELT 

GO TO 149 

00 146 I=1JC 
RAKCIP=RAKUT) FELT 
CONTINUE 


IF (XPS.GT 0.) XPS=XPS* JBL (1) 7363. 


OIA=D1AMN 

ING=IVV 

00 1015C I=1¢ 10 

TLES (192) =INPUT 

TLES (7, 202TLES (A o2)2L INR 


IF( £501) eE0e Oe? TLES (5,2) =LE P32 


IF ( TANS cEQe0e) TLES(6,2)=CALC 
IF UXPS.LTe00) TLES(7 92) =NLINQ 


LF (RAKE oGT De eANDeRAKE oh Eee 1) TLTS (8 _ 7) =NL TNA 


IF(PNOLE* 0.) TLES(992)=PXTNOT 
LF(COcEQe00) FLESC10,2)=CALC 
WRITE (6410066) TOO 


FORMATCLE »4XeI4o* BASIC OF SIGNS) 


LOEN=DEN 
IFCIOENSL T2293 IDEN=4% 
IFC TOENSLTo7) OENSOENS(TNMEN) 


LFECIOENe GE e7 AND UICEQST) NENSOINFIST 


IFC(FHO. EQ. 0.) RRO=1.9995 
IF (SHG EQete)? RHO=149384 
PSC=1644%. 3605 *OTAM*DIAM*2 
PTAN= PCLIZCPIEX3ZCLDD 
AID=N1A 
OUN=9FN 
QAE H=HE ALY 
OHR=REO 
IF(UO.NEeSI) GO TO 15C 
AID=A1LO/ELF 
QUN=OUN/DOSI 
DAE H=QAEH/TLE 
OHR=O0K8/RHOS I 
CONTIBUE 
00 1 [=199 
CAVCI) SVELCD) 
CAE CIP =EHPCT) 
IF (SHP) 103,102,103 
PVEL=VEL (3) 
PCHPZEHP (3) 
CONTINUE 
R(2,15=CD 
BLAS=JPL (1) 
IF (XPS.GEe0.r GO TO 7 
nO 26 T=149JC 
AZZ(T 5 38)=AZZ C1524) 
EXXS=2.*BLAS /PIFSITMPUN (X34 X00 JC) 
FORMAT(2H »2A2,66H 


JEA=0 
IFCEXX(L) NE ele)? GO TO 10004 


W2 


ASKTN FOR" ) 


‘ oat a * 
ES aahmmeiter en tae Tae ee ON 


“ht 


eR ee 


13352 


30 


53 


51 


EXX (LP =EXXS 

JEA=1 

IF(CO.GE.10.) GO TO 390 

TH=B 2-1) 

CFO=.008 

IF (CO.GT.0.) CFO=1. 

IF(CO.GT.-10. .sANO-CO.LT.00) CFO=-CD 

00 10007 T=1,JC 

ITF CCOLES0) TM=10416254AZZ(1 925) 4125.9AZZ( 1,25) 4% 

IF (COLE 710.) CFP=B(T,7) 

B(I,7)=CFO*TH 
DO 15 T=1,JC 
AZZ01,23)=xX3C1) 

00 16 I[=3,JC 

AZZ(Y=-1,19) =X3(T) 
AZZ01519) =X3 (1) 

00 & T=i,gil 

PXTBI (I) =P(T) 
AZZ(1,36) =AZZ(1 425) 
AZZ(1,37) =AZZ(1I,23) 

DO 19071 IE&=1,122 
3(9,2) =JBL (TE) 
XSX=XPS*(360.0/8(992)) 
ASL=XSX/(1e¢0°X3(4)) 
AS2=XSX-AS1 

00 10027) KE=1,IEA 

TLES€3,2) =TLES (4,2) =INPUT 

EaPU=EARF 

ERA=ABS( EXXS-EXX(KE) ) 

I= CERALT 2.005 .AND. KE. GT oe Le AND. JEAEQe3) GO TO 10070 

IF CERA.LE.2005) EARU=EARD 
00 19069 IRP=1,IRPM 
KI=0 
RP M=XMM(TRP) 

EAR=EXX (KE) 

00 100 I=1,J/C 

3(1,3)=X3 CT) 

91194) =X4(T) 
B(I,5)=x5(T) 

O00 19051 LE=1L, NC 

B(LE,6) =(BLAS*E AR*¥X6(LE)) /(B (9,2) FEXXS) 
3DEL=A9SS( B(3,6)-X6(3) ) 

IF (BDEL.GT.»c0001) TLES(3,2) =ADDU 
00 19252 LE=1,JC 
AZZ(LE,25) =AZZ( LE, 36) *B( LE, 6) 
JO Sf T=159 
VELCI)=CAVIT) 

EHP( I) =CAE(CT) 
8(5,2) =XMM(IRP) /66.0 
IV=i 

IF (SH®.NE. 0.) GO TO 51 
00 53 I[G=1,IVV 
VELL CIG) =VEL(IG) 

EHP4 (IG) =EHP(IG) 

GO TO 214 
VELI CL) =PVEL 
EHP1(4)=PEHP 


ara) 


21 


116 


15 


10035 
101 


197313 
10719 


34 
10045 


N0 5 I=isit 
AZZC1, 249 EAS1 *X3 01) 4AS2 
SKZ=ABSCAZZ(1.26)) 
IF(SKZeLT2 0001) AZZ(1,26)20. 
IFCXPSeLT: Ood AZZITe2od=AZZ0 Fe 38) 
AZZUL, 38) =EAZ7 (1425) 
P7701, 25) =AZZ (1, 36) 
A7ZCL_ 230=A2Z (1937) 
IF (SHF.E0s0.-OReIVeEQei? GO TO 18 
NO 114 JeieIVV 
8115) =VELCS) 
A3 (J) =FHP CS) 
S1=VELACIV) 
CALLDISCOT (S19 S198119%3563 9-120 550 00yS2) 
THPL(IV)=S2 
B(Ey2) 2 (325.86%EHPLCIVID/ (VELL CIV) SETHEUS) 
P(7,2) 21666789 VELICIV) ° 
NN=6 


0OJ=1.¢ 
JEESTANSI 
J00=004 
JO=- €666E7*F_ OAT (CJC) 
R(4,1)=2.6¢ 
R(2,1)=CO 
21351) =TANB 
B(6,1)9=0.85 
R(6_1) 2920 
P(7,1)=20.0 
4(8,1) 20.6 
P1(9,19=C.6 
P(4,22 21-0 
6(2,2)=1.0 
9(@,2) =P HO 
A(Ssl P=EWAKE 
A(3,2)=05A 
(4,2) =H" AD 
PSL=B(7 4297030 14159265 9B 15,294 2(3y2)) 
LF(TANS.GTo0-) GO TO 101 
NO 12935 I=1,JC 
&(L,2) =2SL/B(1 93) * 80164) 
AJJ=12C€ 
IF( 3(195).GToce) GO TO 10019 
IFC BUSolP clEeDe) BISolL VHA II04) 
00 16018 IsisJC 
BCT, 5)= RSLESQRT (B15 etP FBT eGVIS (SUT, 32 F304 01)) 
NN=NN#1 
PAKOP=L 
LE CRAKE LE ee Otc AND RAKES GTO06? GO TO 190% 
00 34 I[=1.i1 
RAK (TPS RAKE (XS CTPAKSCLIIS (KS CALF XT (1)) 
VK=8(7,2)/1.6878 
IF (KIeGT.0) GO TO 85 
AJJ=16€ 
WRITE (6510025) 
PHS=SHP 


a tenes ena 


re 


Decne nuchuaenwotcme stasis 80, etery eT 


00 152 I[=1,Ivs 
IF (U0.-EQ.SI) GO TO 151 
OMVCIDEVEL (I) 
OME CIYZEMP IT) 
GO TO 152 
4151 OMVIII=VEL(LI/VST 
OME CI} =EHP (13 /PWR 
152 CONTINUE 
WRITE (6.13008) UI,U0 
TFCSHP.NE0.) WRITE (6,10021) PU, PHS) ,OENL»CENZ.NEN3,0UN 
IF (SHP. E00.) WRITE(6,10020) SHPT1,.SHPT?,DEN1,.9°N2,G2N3-DUN 
WRITE (6410027) VL eVU, (COMVIT) os I=1,I1V9) 
WRITE (6019045) PE,PUs (OME (IT), Tz1,IV9) 
WRITE (6210028) ULe AIDE WAK= SC THRUS, UL NAZEH,UIFN,OHR 
19525 FORMAT (72X%— 90% AG, F= Se FB HDF ato WTT=% F564 4F o1-THD=4,FS5.4,*% gH, 
2 Rag t=% oF Boee® eRHO®, A110, FEF, 1X oF G0G5/) 
WRITE (6,1903C) (JBL(I),I=1,I27) 
WRITE (6.10031) EXA, CEXXCID,T=1,IFA) 
WRITE (6,10032) URPM,CXMN(I) -T=1,IRPM) 
10025 FORMAT (141) 
10523 FORMAT(/,287,2A10,A8eF10 04,7) 
100214 FORMAT(/ 9? PS (4% A359 FE%_ 1PEL2. Uy 2ALG AR, IPFINWG LS) 
19633 FORMATC2X,*2%,9X,9I3) 
1DO0T%1 FORMAT (2X ,AEe SK, IPIEL2Z.%) 
10032 FORMAT(2X,*N® ,A9,1P9E12.4) 
10627 FORMAT(2X%,42,A6,2X,1P9E12.4) 
19745 FORMAT (2X%_2A344X,_,1P9E12.4) 
WRITE (6,10200) (TLES(I,1) .I=1,10) 
WRITE (6010161) CTLES(IT,2),T=1,16) 
10103 FORMAT(/ 1X, 1003X,AL0)3 
10101 FOPMAT(C1X, 1003X%,A10)./) 
OO 19111 T=1,J/C 
IFCFOO.LE.0. PITI=0 
SRKZRAK(TI/IIA/S12.6 
19411 WRITE (6410102) K3ICT) ge XS CTV S01 96) AZZ 61425) SUT eS de BUT 94D 
2 AZZ(1 424) sS2KyP (I) BCI, 7) 
13162 FORMAT (1P10=13.4) 
485 90 1€090 [=1,JC 
A27(1.25) =AZ7 (1,25) *BCT, 6) 
13095 RCI e3IIAB(I 25) FAIS 
IF(JEE.LE-0) GO TO 4&7? 
00 209 I0C=1,3 
NO 201 T=1,57 
21 RIIT,30=3(1 30) * ASS 
9146(61)=.975 
814(62)=1.000 
814(43)=1.025 
NO 215 [J=41,43 
TJTSi Js +& 
TUP=IU+8 
00 216 T=1,.JC 
216 Q(1Te5)29(7, 30) #81401) 
CALL SUB 
ALGCT JT) =PP7 
3144 (17 JP) =PPB 
2145 CONTINUE 
JK=64 


q TTT=B(6929/1(B (89 2)*B 0392) 982930141592 79R 1747) %92)/R) 
i] IF(€ Bliy2)05T.0. } GO TO 16 
% JK=&GB 
TIT=TTT*550./BE7 02) 
10 M0 11 T=193 
CCl, 1) =BiGC KET) 
00 11 J=103 
K=3*(J-1) I 
B15 (GL#K) =BLb (UC H) OF (T-1) 
11 C (Je 1) =815 (Kt60) 
CALL SIMEQ(3+0,CC) 
20s AJSI= (“CO (2921) FSOPT (COE 2g LIFE PHU FEC (Fe LIF (CTL eo LIMTTTIIIS(2.9CO( 3, 
11)? 
48 Wy G9 1=1,JC 
69 95195) =8 01,30) FAIS 
G7 JEC=0 
CALL SUA 
PP11=315 (181) SET HRUS 
PE12=cHPACIV)/PPLL 
THP=2 (6,42) 
ATHR=3(B 92) /302¥ 30141598 (R(Z,2058 294 (R (7, 2) 842) 2 OPT 
ASHP (IV) =PP1i2 
TFCSHP) 55910950 955 
g 55 CONTINUE 
KI=KI#1 
IF (K1.GE.63 GO TO 10050 
IF (ABS ((SHP “ASMP ITV) D/ASHP (IVD 969505) 12959619055 42¢ 
7 2) NI=.33 
{ IFGIV.EQ.1) GO TO 822 
NIZALOG( VELACTV<LD/VELE CIV) }ZALOGE ACHOCTV=4)/ASHOCTY) ? 
823 VELICIV#LD=VELLCIVD*C (SHESASHPCIVID SENT ) 
AVLEVELACIVEL 
LE CAVE GTeCAV(5) oORSAVL oLTeCAV(12$ WRITE(E,927) AVL 
922 FORMAT (////19%_*ESTINATED VELOCITY VALUF IN IT=RATION FOP OFSIPED 
1SH® IS NOT WITHIN RANGE OF INPUT VELOCITY VALUES. 0e®/15X 0% ee ePROGR 
2AM CANNOT EXTRAPOLATE FOR COROESPONDING FHP, ESTIMATED VELOCITY 
3VALU: = *1PE10.4) 
TV=IVed 
GO TO 21 
12250 0 12669 IX=1,IV 
FHP CIX)SEHPACIX) 
100469 VELCIX)SVELL(IX) 
00 4CG50 I=1,C 
AZZ(1 926) 281914) 
AZZ(19279=3(195) 
BI=O22(1y27) 
IF(P00.GTe 00? BI= PUID/PI/X36T) 
CIISAZZ (1924) ®B(Ie3)*SO9T (104 TI*BI) 297/130. 
AZ7(1y2802CII-S(T 96) 
AZZ(1,29)SCII#¢B(1 46) 
AZZ(1,3C)=B(1 96) 
AZZ1y 21028 (19 4) 
AZ7(1932)=B (1912) 
AZZ(1, 33028 (L913) 
AZZ(1 9 34) =AZZ (1, 250% 200 
49950 AX(T)=AZZ(1934) 
a. IF (AZZ (11925) eNEo Qe) GO TO 5555 
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& 


& 


SUP CAZZUD) BH) AZZ (65 34) ) JAZZ 095 23) AZ (622 99) 
YINT=AZZ (99 36) -CAZZU9e 23) FSLP) 
AZZC10 BUD =ZCSLPFAZZULO 230 PFVINT 
AZZCLLs BHD ECSLPPAZZULLe2SPIOVINT 
5555 DO &6952 K=26, 3% 
00 40051 J*i,%1 
BLi105) =B(5, 3) 
0051 B3CJPZAZTIS eK) 
00 43052 I=i,11 
S1=AZ7Z (1223) 
CALL DISCOT(S1 »S1 0811983 993—-420 9 IC 0952) 
0052 AZZ(1.KI=S2 
00 SC5C T=1.JdC 
AZZ(1,2) 28192) 
POT, 291=H801915) 
Hil ,2L)=B01 05) 
BI=B(Ie2i) 
JF ¢Pp00.GT.c.) BI= PCLISPISK3(T) 
8(1.22)=8(1,38) 
B(1,.23)28 (105) 
R(LTy24)=AZZ01, 28) 
5([e25)2A2Z2 1925) 
TX(TI=SBUT»25) SOTAF1200 
505) B(1e2) HBT, 24) * BCT o3b*SORT (1. BIPBI) O11 95~ 
AT (11) =B01145) 
OPRTAIN DATA AT 2 , 
00 £0680 [=1,10 ae 
O=AZ7Z (1,19) ‘ 
GO $0239 K=20-24% 
0030 CALL PTSCOT(D yD oBlLe3) BULK) BCL eK) g— 1799 Io De AZZIT KD) 
CALL OTSCOT(D sD eBlin3) oP ePe 120 yJ0o0 CX)? 
CALL OLSCOT(D Do Blt s3)eBlTe18) BUT g 18) or 120 IC oo AZZCT 9180) 
0:82 CONTINUE 
Do £0090 f=1,10 
07935 P&TP=CXET) 


CALL STRFSS(EZZ,AREA eXBA Dy YBARs, AVEXO g AYEY Dg SMX pF MYO yEMT Fg EMN% STO 
1MAX RAK, PAKOP) 


C THE FOLLORING STATEMENTS HAVE REEN ANMED IN ORDER TO SEND THE R 
C VALUES TO SUBROUTINE WEIGHT. CHORD e THICKNESS + AND CAMBER ARE fT 
Cc PITCH ANDO SKEWR APE IN RADIANS» 
VS=VK*1.6878 
nO 999 T=19J0 
IF (PD0.GT«0.) PUTY=PXTBICT) 
KTVT. 3) 
IF( XL) eGTe 0665 eANDe KCL) el Te 075) SIGMASR(I,19) 
CHOPD( I) =BiIy6)*O 1AM 
THICKNS (I) =AZZ (1, 36) *DTAM/2. 
PITCHCI) =ATANC BIT 5)? 
BBUCIPZATAN(3 (Ty 8)) 
FX(TI=ACSG 
IF (CLIC 1) eE2e Of) GO TO 29 
EX (Late O/ Cte OF (Go 28329TANIPITCHETI-BRICTI) SOLA ITD) 
29 CONTINUE m 
QIEP CTD /PI/K SUD) \ 
LESPOO.GT.0.) PLTCH(I)ZATANCOT) ‘ 
SKEWRS TY ZA7ZUT 9 38975722958 
117 
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OZZ7(1 SHV =AX(T) 

AV =(39(I,6)¢B(1,129)9%2 

BV =(9(Te4)/B(T o8)-B UT. 1 329 9F2 
VSUBFE SO CTI 2VS**2° (AV *8V ’ 
VSUBRCT) =SORT CVSUBRSOCT)) 

999 CAMBS RIT)=.06799R (1,18) SOT AM 
PPS=RP M60. 
00 996 Iti, 
PXTSICTI=PL83(1, 309861 5) 
PXTB(IP=PIFBI(I 43) SBC, 83 
936 CONTINUE 

EALZ CE ARTI. 1G1SIFOTAPF2) 600 
AL=3.16159%0. 7°B (7,5) 
AP=EA1°(1.967-P.229°AL) 
VA=VSFRI7,4) 

VR=SORT(VAF92¢ (0.793.161 E99RPS POI 1D FF 2) 
TC=220%8 (692) 7 (RHOF APF VRE 2) 
SIGMB7= (64.4° HEAD) /VRE°2 

WRITE (6510242) 

10242 FOPMAT(1H ) 

UVE=VU 

IF(CUO.NE SSI) UVREUVE 

WPITS (6,78) UVP 

90 75 J=1e11 

VPST=VSUERI{T) 


IF(JO.£9.ST} 
79 WRIT? (e860) 


VRS ITzVRSI%. 306A 
BCL ode BCL Se BIT Bde BEI old yg FIT of Pde FCT e120 69(T 1%) 


1 sBCTo179eVASTe30T,19) 


WRITS (6.10042) 
WRIT: (6,77) 
WRIT! (6561) 
WIT? (6410042) 
ULO=ULE 
IF(U0.£0.S7) 


PP1,PP3,PP&,PP2,PPS5 
PPE »PPBePP9,PPT,PP1) 


ULO=ULS 


WRITE (6476) ULO 


MO 76 TH1sit 


SLI=.e*RHOFVSUBRSOCTI*CLI(TI® 81.6) "NFA 


IF (UO.EQ.ST) 


SLI=SLIF14.5939 


75 WEITE (6520065) VCI1e3) -CCONE CTY, COTWOC TI» TCTHR(TI Co SGOUITD oFX(TD 


2 SLI yAZZ(1 438) 
WRITS (6, 10942) 
ANY="AV 
Pi2=PP12 
RTH=THR 
UB=UFT 
UTO=UFE 
BTHF=ATHR 
UST=UPIN2 
VLS=3 (7,42) 
V2=UVE 
U02=UL¢c2 
YOG=ZULZES 
UMO=UIL 
UST®=UL PIN? 
IFCUO.NE 6ST) 
uN2=ULS2 
UO4=ULS& 


sAZZ( 1928) pAZZIT 929) o AT7IT, 34) 


GO TO 16C 


UMO=UNM - 
USTP=UPA 
UST=UPA 
UBSUSL 
UTO=UFS 
V2=USV 
VES=VLS*.3068 
P12=P1i2/PHR 
RTH=RTH/URT 
BIHR=8THR/UST 
160 WRITE (6,10044) PP11,PU,P1L2,.ETHRUS, EWAKFE, VL »pUEV, VK UTO RTH, USL (IE) 
2 eRPMHEARUSE KASEAR VL oe V20 Vi SSUTO-STHR 
TFCUO.EQ.ST) HRITE(6,709019! UWL2 UBeUB 1104 104,UM0,UMO0,UMO,UNO,UST 
IFCUO.NEeSI) WRITE (6,200093 
IFQCUOLME.SY? HRITEI6,299263 YO? UB, UE e 1iI0&,U04,U40,UM0,UM0,U40,UST 
77?) FORMAT (GX, FOP TY=%_ 1 PE10 eee FCPSI EF LPEL0 e Uy GX oe FETAL HS IPELC oly 
2 &UXyFCTSTE* IP F102 GK sFCTSIEZCPS IF gi PE10.&) 
75 FORHATCOX,* MF eg OK gFTANBIF 7X 5STAN BE 5X _ FG4, 9K FUT/S2V4_ 7X, 
2  FUASZSEE 7 Xe FOC TSI 4g 7Xy FOCPSTe pS Xe PVR IF GAT 5X, FCAVVE) 
639 FORMAT(1P10E12.&) 
81 FORYMATISKX,*#CPT=*,iPEL0. &, SX p*OPSHh ey LPELG eb y5Ky *ETARH%, IPELC 1G y5X,y 
2 *CTS=*,1PEIC. Ge OKs *CTS/CPSH* 1 PEC 4) 
75 FORMATC5Ky* X¥s8XoFCL MeO FALT (DEG) FISK, FFM/C#. 7X ey FCDO/SCLY 7X, 
2 FE CXP SS aXe ALITY, AB 2K, FTETSEGEGE FON CF (C/PODILE* UX FICS OD) TEP EX, 
3 *T/F0F } 
25546 FORBAT Ui PEL3.3,2P10541.33 
10544 FORMAT(E1X%, FETAD=F,1PE10 oe Ho GX a FPS (Fg AG, P= Fe LOPE LS og 3K FL -THN=*%, 
ei PE10.4% yg pas 
SX eF 1“ WTT=H*_ LPE 10s by 3X AZ, AB yg *E*o1PE10 Gy Xp *DOSIGN TH*,AS,%=%, 
LPL LG eG eg GK gt Z=FT2— BXyFNCREVIMINIRERF,1PH10 Gel EX, 
AB AS, F=F,LPEL0 ee SKeA2,A7, * =* gf PE1964.3X%,*CALCULATED TH*, 
45 ,*=*,1PE10.4%) 
10 746 T=197 
SX(1?}=4ZZ(1,19) 
APR=ARVEA(T) 
XPP=XBARQI(T) 
YPR=YRAR(T) 
XOT=AYEXOCT) 
VYOI=AYEYVO(T) 
XOM=EMXO(T) 
YOM=EMYO(T) 
TAM=EMTB(T) 
QBM=EMG8(T) 
STRM=STQMAX(T) 
TFCUO.NE.SI) 60 TO 74 
APR=APRFEL2 
XPR=XPR/ELT 
YPP=YPR/ELI 
XOI=XOI*EL 2° EL? 
YOI=YOI *EL2*EL2 
XOM=XOM/SIM 
YOM=YOM/SIM 
TEM=TRM/SIM 
OBM=N8M/SIM 
STRM=STRM*68 94.757 ~ 
74% WRITE (6,20047) SX(I) sAPReXPRAYPR»XOlI se YONI, XOM,y YOM, TAM,QBM, STR 
20067 FORMATCLPE10. 391P5E11.3,1P5512,3) 


Mewar 
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20009 FORMAT CIHL GX oe FX 4, 9K, FAREAF, 7K yg FXBARS 7K, FV FAQS 7X, FIKO® AX, *F IYO” 
2 0 We FMXO% 4 9X o FMYOF 2 9X, FMT AF LIX, PMOBE BK, PMAXSTRESS® 
ZONLN FORMAT (15% 05 OX pAb 7X pA bX yp AG eg BK AS AKA F UK yg AB GK AB UX AR GSK 
2A9) 
10019 FORMATCLHL GX op FX Fy 7K y FAREAS AS UX eg FXAAOF AG AK FVGAPH,AGS2K, FIX 
2  —_g@hd oe IK eo FT VO™ A5 GX a FMKO% AB a 1X yg PMVOF ARG LK FTA ARG 1X, MOSS AB, 
3 *MAXSTRESS*®,A9) 
WRITE (6510042) 
WRIT: (6516042) 
WPIT: (€,106011) 
00 73 T=i1,.i11 
SRK=RAK (I) “JITAS12.~ 
WRITE (6519012) BCI eS) SRK PXTSICTII,PXTA(T) 
73 IFCPOO.GT.06.) PXTBICTI=PCT) 
19621 FORMAT(5X_*X¥,8X,*RAKG/DN® 1X, 23H PY YTPANIIT PY YTANA ) 
10512 FORMAT (1PE10.3,1P3£11.3) 
CALL WEIGHT (JC eS IGMA7,HUB,PMHC,. WEIGHTR WE TGHTH, FAK) 
SQI=SINC PITCH(1)? ?) 
321% HUBSPAC=2.079T FX (1) SOT AF6.90/2(952) 
TAI=AY(1) "OL AM*6./SBI 
RL ASPAC=HUBSP AC-TSI 
FILS PACH=BL ASPAC- OF TRI 
WRITE (6,998) TC,SIGMA7 
BLOE=38LASPAC/DIA/12. 
FILS=FILSPAC/OIA/IL2. 
UF=:FT 
UMI=PMI 
IF(UO.NE.SI) GO TO 168 
UF=U0Q«6 
UMT=Uc7 
153 CONTINUE 
WRITE (6,994) BLSe FItS 
098 FOFMOST(/20Xs*S3URRILL THRUST COTFF TOES TLV G//20KX,*#8U9FTLL CAVI 
LTTATICN COEFF SIGMA(0.7)=*,=1C.4) 
994 FORMAT(/20X,*CLEARANCE AT HUB RETWEZIN SLANTS/SI=* LOE L 2G eS e 
2 2C¥,*CLEARANCE AT HUB RETWEEN FILLETS/92¥%,1°F 12.4) 
nO 122 T=1,l11 
ARPMT (T)=39 (1 96) *AZZ(1,25) *PSC 
122 PMOXC(I)I=ARPMI CID FXO) FF2 
PMOFIA=9 (992) SDEN/S1728 6% (OIAMFS OD EFSESTMOUNIK,OMOXK, JC) 
TP MIN=PMOF IA P MHC 
PADOGB=SORT (PMOF I B/KREIGHTAI/12. 
PANOGH=SQRT (PHHC/WE IGHTH) /12. 
FANOGT=SORT(TPMIN/ (WE IGHTBtWETGHTH)I/12. 
PMS=PMOFI8 
TPS=TFMIN 
GRB=RADOGS/ITIA 
GPH=QRANOGH/OTA 
CreT=RADOGT/OIA 
IF(UO.NE-SI) GO TO 171 
PMS=PMS*.00029264 
TPS=TPS*.000292E4 
171 WRITE (6,995) UMIT»oPMS,UMT, TPS,GP1,GR2, GP3eGWL +622, GRH, 
2 GR1yGR2, GRT 
9395 FORMAT(1H1e19X%_*MASS POLAR MOMENT OF INEOTIA OF BLANES *,41), 
2 €F13.697/7,20K% ¢* TOTAL MASS POLAS MOMENT OF INERTIA F, ALC ELI be SS, 


3 20%,2A9,* OF BLADE/Oe* FG eS/20X%0249,* OF HUS/D2* aE Veh, 
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uw 
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wn 
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$95 


596 
597 


519 


& 4/,20X_e*TOTAL %,2A9, /02*,F9.4) 
HIS SECTION CALCULATES THE ABS COEFFICIENTS 
IF(TYPE.EQ.6.) GO TO 579 
PAD=T.25 
IF (8 (1,3) .6E.0.25) RAO2B (1,3) 4.05 
IF (8(2.3)-EQ.0.25) GO TO 599 
00 510 J=1ell 
B11(J5)=R8 (5, 3) 
SI(JP=PXTBI (CU) 
$1=.25 
IF (8(1,3)-6E.0.25) SizRAD 
CALL OISCOT(S1,S1,8114,83¢B3_~1204 JC .04S2) 
PTWOFIV=S2 
0O 513 J=1,11 
§81205)2=8(5,3) 
B3(J)=8( 5,6) 90 ITAMN*%12.0 
$1=.25 
IF (8(1,3).GE.0.25} Si=RAD 
CALL DISCOT(S1eS1,611,»5 834.83, -120, 50095S2) 
DJURL U=S2 
N0 5165 J=igii 
BLL JPA, SI 
BIC JI SARPMI (5) 
S$i1=.25 
IF (8(1,3).G6E.0225) Si=RAD 
CALL OISCOT(S2,S1 4811983 4834-1205 IC 405S7? 
ATWOFIV=S2 
00 5€9 J=1,il 
Rii(J)=3 (05,3) 
R3 (J) =THICKNS (5) °12.0 
S$i=.25 
TF (6(1,3)9.GE.0.25%3 SizeRAD 
CaiiL TISCOT{S1 $1,811,483, 83,-120 55049587) 
TMAX=S2 
UF=6.5*T MAX 
GO TG $92 
PTWOFIV=PXTBI (2) 
OUBLU=8(2,6)* D0 ITAM*%12.0 
ATWOFIV=ARPMI (2) 
THEX=TRICKNS(2)*°12.0 
UF=C.S5*TMAX 
CONTINUE 
IF (8(7,3).€0.0.7) GO TO 596 
00 595 J=ieil 
§11¢5)=8 05,3) 
RI CJI=PKTBI (J) 
Si=.7 
CALL OISCOT (51,51 ,811,83,83,-1204 JC e606 5S7) 
PSEVEN=S2 
GO TO $97 
PSEVEN=PXTBI(7) 
CONTINUE 
RAKE 2Z=PAK(11) *DIAMFAZZ (12,38) FP (1LI/3C. 
00 519 J=1,7 
B1idUP=SXCd) 
RICJPSAVEKO (J) 
$i=.25 


oA 


oe 


$35 


C43 


572 


IF (SX (1) .GE 20-25) S1=SX(1)40.95 
CALL DISCOT(S1L .S1.B11,83,93,-120,7+0.82) 
AYO=S2 
F=INT (TYPE /2) 268.4 INT (TYPES3) FS. -INT CFYPSS4) 956. °INT(TYPE/6)%127. 
WL=INT (TYPE/ 2) * .SHINT (TYPES 3D FC OLS INT CTYOE/S ODF, S2-INT (TYPS/6) 230 
CS=ATWOFIVZ (OUBLUSTMAX) 
CN=AYOZ (UF FJUSLUFTMAX®*2) 
AA=1.0#6.0/PSEVEN+G. 3#PTWOFTY 
BA=4300. *RL FEAR 1B (95.29% (RPOM/100.) FPF (FIAM/ IDL) FF 
CA=(1.041.5*OTWOF IVI * ( OUBL UF FHA) 
ANEWE13.*SO2T (AAFPPL2/7(°(CAYRPMFR(9,2)) ) 
BNEW=CS*BA/G./CA 
BTHICKL= ANEW/SORT (CN) + BNEW/CNFRAK (11) 
RTHICK2=ANEWZ/ SOR TICN) CBNEW/CNY RAKE? 
WRITE (E4530) RAD 
FORMAT (////29X% ¢* ABS COEFFICIFNTS (CALCULATE) AT THE*,F4.2,* 2ADIUS 
19") 
TABS (£2=TBS 
TASS (6) =UBS 
TFCPOO.GT.0.) TASS (5)=VBS 
BSL=ATHICKI/OIA/12. 
BS2=8THICK2/DIA/12. 
ATF=ATWOFIV 
IFCUO-NE.SI? GO TO 172 
ATFSATF* 00064516 
WRITE (6,561) TAPS ,BS1»BS2eAh sg BAgCAGCS oONySNIR,ATE 
FGPMAT(//20% ,6A10,% 


2 /cCOX*USING ABS RAKE = CONVENTIONAL PAKE, T/D= *E10.4/40X*UST 
BNG ABS RAKE = CONVENTIONAL + SKEWeINDUCED PAKE, T/I= FELC& 
2 //22K*VALUES USEN IN NETERMINING THICKNESS= 
3 A= *E12.6/62X*8= FE12.6/62K4C= FE17,5//°CX*SECTION AREA COEFF 
GICIENT C3= *EL1064//20X*SFSCTION MODULUS CIFFFICI-NT CN= F254 
SleU//2CK*AREA OF EXPANDED CYLINDRICAL SSCTION IN SQe*oh695X_FAS=F, 
6 £10.4) 

CN=,1 


BTHICKi= ANE W/SORT(CN) + BNEW /CNFRAK (11) 
ATHICK2=ANEW/ SQPT (CN) BNE W/ONFRE XE 2 
BS2=BTHICK2/0IA/12. 
BSL=BTHICKI/DIA/12. 
WITS (6,543) TABS s8S1,AS2 
FORMATIS/S///7420Xe*FOR CN=e1*o/ 20K GAIL, 
2 /4SX*USING ABS RAKE = CONVENTIONAL RAKT, T/M= FELGC4/4OX*UST 
BNG ACS RAKE = CONVENTIONAL + SKEW=INOUCED RAKE, T/N= FE1C.4) 
CONT INUE 
BLADE L=A8 (1153) -B8(1,35) 
CEL TA=PLADEL/S. 
NELTAL=DELTA/2. 
00 i132 [=1,11 
XR(T)=B(153) 
TANBETICI)D=6(1 45) 
GCTI=B8 05,14) 
XSLCTPF4Z2Z(01, 25) 
XSTCIVSAZZ(1,29) 
WAKE (CTD=B (194) 
UVELA(I)=B8(1,12) 
THICKR(T)2AZZ2 (1934) 


eT 


132 


130 


131 


133 
1 


13% 
2 


174 
135 
1 


175 
136 


100659 
10070 
10571 
10672 


CONTINUE 

KR(ILZ2I=KXRILIOCDEL TAL 

XRILIVEXREILI*DELTA 

00 130 I=14,19 

MROTIPEXR CISL) DELTA 

CONTINUE 

XP (26) =XR (19) *OEL TAL 

00 131 [=12,20 

CALL DISCOT (XRT) XRT) go XR oT ANBETI 4X2 4 -044511,0,TANSETICOI)) 
CALL DISCOTIKR EIT) oXRIT) XR oGoKRy -OGGe 1150eG6(T)) 

CALL DISCOT(KR IID XRT gp XReXSLeXRe-Ohbet1eD XSLT)? 

CALL DISTOT(KRIT) XRT) oe XRoXST eo XP e-Ob ell yO eXSTLid?) 

CALL DISCOTEXRIT) eXRUTD a XR gpWAKE 5 XRy Dg tig J eWAKE (ID) 
CALL DISCOTEKRIT) oXRULTD gXReUVELAGXRe 0451190, UVOLACI)?) 
CALL DISTCOT (XAT) pXRIIT Ve XReTHICKRe XR, “04601 Le0gTHICKRII)D) 
XSL CID =XSLCI) POTA%6. 

XSTCI)D =XSTCI) FOIA%6. 

THICKR CT I=ZTHICKROIID®OITA%6, 

CONTINUE 

WRITE (6,133) 

FOPMAT (LH1—9////e10Xe*RADTAL PROPELLER DATA FOR INPUT NTO DESIGN P 
ROGRAMS(§ RADIAL STRIPS ASSUMED) F,S/7) 

WRITE (65134) SNIP,SNIP,SNIP 

FORMAT(ICKs® KR¥gSXe*TAN BETA [99 BX oFG4e7K oe F*KSL (4 A7 SX FXST (GAT, 
SK oF LOW, BK 5 FUASZVSS SK yg *THICKNESS(% Aly /) 

00 174 [=12,20 

SSL=EXSLCT) 

EST=XST C(I} 

STHR=THICKR(T} 

IFCUO.NEeSI) GO TO 174 

SSL=SSL/ELTI 

SCET=SST/ELT 

STHF=STHR/ELTI 

WPITE (O9135) XRITISTANBETICIY, GEI) sSSLySSTy WAKE (ID LUVELA(I) »STHP 
FORMAT (GX F106 o So SX oF 10 So SKE Ze Gg SK FLV Fy GK gH ee LK oF LC eS gy bXyel 
fey SK,F 10-5) 

XSLCLLI=REXSLILL IF OIA*E. 

XSTCLLINAXST ULL IFOIASS. 

SSLEXSLLL) 

SST=XSTCL1L) 

IF(UO.NEeSI! GO TO 175 

SSLESSL/ELT 

SST=SST/ELTI 

WOITE (6,136) XRC11)-SSL,SST 
FORMAT (4X —F10059 SIX F100 59 4XeF 106.5) 

IF(CSHP.NE.0.) GO TO 10069 

IV=IVed 

IFCIV.LE.IVV) GO TO 2t 

CONTINUE 

CONTINUE 

CONTINUE 
CONTINUE 
STOP 
ENO 
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SUBRGUTINE SUB 
OIMENSION CHIVD(LL) »>THICKNS(11) ,CAHMBERI11),°TTOH(11) SMEWROLL) 
1 eXI¢11) 
OIM~NSION P(13),8T7 (11) 
COMHON/C WEIGHT /XI » CHORD, THICKNS » CAMBE ©, OTT OH, SKEWR NIA M77 50EN 
1sPAKE sPI 4g PP7, PPB, PPG, PP11 gy EWAK® pVSeRPS SIGMA, EAM, BT oP, POO 
2 oF WOCTAM,AT TOI AM,HUBLEN»FORORE ,ADBORE,NISQF! 
OIMENSION 8(38,38), B2C1 RAL) eB3ICL481) 
1 ¢B1i0181),916 (181) 
2 0GB15 (181) pAZCLLsgI1L) pBR(L Le itd pCOL2 263 ,50C(12,12) 
COMMON 8, B2.B3>5 Ril 981 4,915 ,A7,5H, 
1C,CC, 10s J8eJC JO ,JOOSIES oe CLIC11) 
COMMON CCONE(11) ,CCTWO(11) ~CCTHO(11),CCFOR (11) 
COMMON 921,P92,PP3,PP4,PPS.PP6,PPIC 
COMMON SNO73),CS(73}3 
10252 0 10625 N=1,JCe 
NO 2C021 T=1.JC 
ABG=1./9(1.5) 
AAH=B(Ne3)/B(T,3) AAG 
AAQ=B (1,5) 
IF CABH~84G6)16619, 10018216919 
1321a B2(1)21./SORT (1. +AAQ*AAQ) 
B3(1) =AA0*B2( 7) 
60 TO 20021 
109219 S=1.+AAHe%2 
T=SQ=7T (S) 
V=i. tAdGs*%2 
W=SQ2T (Vv) 
NE=T-wW 
U=EXP EHO) 
R= (00TH1.)/A8H® (AAGS (WH, 20) PU) 991942) 
80=.25 
X= (107 620%8(9 92) SAAGIIFEC (VIS) FFA0) 
Y=(09, FAAGEF2D 2d S(VFFACD FCI S.SAAH FF2—-2,)/1S*FAC)) 
7=1./0246.9309,2)5 *Y 
IF CAAH-AAG) 1990 21,10021,10023 
10820 AFHL1.+1./(P-1.) 
AA=X*(1.70R-1.)-Z5AL OG (AFD) 
P2IT) =2e°B(9,2) F295 AAG AAHT (1,-AAG/AAMD FIG 
AZ(I) =909—92) * (1 ee AAG /SAAH)4 (1,426.99 (9,7)"°AAG FAA) 
50 TO 2002L 
1D52L AGL. Fi eSUL SQL.) 
OB=eKF(LSILe/WLedDHeZPALOG (AG) ) 
A2(T) == 9.0992) SAAGT (1. -AAH/AAGI 4 (1,22. 499(9.2)* AAGFAN) 
331) 8 =2.9°B19,2) ##29AAG* (1, -AAG/AAH) #A8 
20021 CONTINUE 
20524 FOPMAT (9F1206) 
00 2 1=1,JC 
2 72107) 282,39) 
no 3 [=1,.37 
$1 TeF5F(16 4961, DHF 301 SIIFCS ET) 
CALL OISCOT(S1,51,811,863 33 4-120¢3040,S7) 
3 R1iS(1TI=S3 
N0 5 [21,37 
Si =o57F(1.tB 01 SPM EFC. 9C1e 3d) FCSCT) 
CALL DISCOT(S12sSi»P14982 632 92129950 eo fy S) 
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5 B14(1)=S2 
00 & [=1437 
B2(1T)B14(T) 
& RILID=B15 (1) 
00 10022 L21,35 
NL=374+L 
N2=37-L 
B2 (NL) =B2(N2) 
10022 B83 (NL) =83(N2) 
C2=2./726 
NP=72 
NH=36 
XNP=NP 
S=C0.6 
SL=0.0 
DO 20 I=1_NP 
S=S#¢B2(T) 
2) SL=SL+B3(T) 
Bt(1,9)=S/XNP 
R(1,10) =SL/XNP 
N0 &¢ T=1,_,NH 
S$=0.6 
St=G.0 
BO 32 J=1,NP 
K=(J-1)*I 
K=MCD (Ky 72) 61 
S=$+#3205) *CS CK) 
35 SLESL#BICIPFCS(K) 
L=T¢+1 
Bt(L,9)=S*C2 
43 BL,1C)=SL*C2 
CPHT=( (1 e438 (1,3) -2eFBING SIPS (L.-BU1,3)) 
IF(CPHI «LT eis) CPHI=-1L. 
LF(CPHI.GT el.) CPHI=1. 
3(N,11)=ACOS(CPHTI) 
RC1,11)=.-6 
B( JC 11) =3.14615927 
CON3=3.1415927 
NO 10025 T=1, IC 
SMP=SINCFLOAT (CI) *O(Ng11)) 
CMHP=COS(FLOAT (IT) *B(Ng114)) 
IF (N-1910027, 10026,10027 
13027 ITF (N-JC)10028,10€29, 10028 
10926 AZN=.6 
AZN=.0 
N2=I +1 
00 20626 K=1,N2 
IF (K-JC) 10070510070, 20026 
1007] AZN=AZNtCONS*FLOAT(II*BIK, 9) 
BZN=BZNtCONS*FLOAT (ID *B(K,10) 
29026 CONTINUE 
AZ7L=.0 
BZL=.0 
IF (N2-JC310060 ,10030,16030 
10662 Ni=N2¢e1 
NO 20036 M=N1ieJC 
L=M-1 


ety} 


20036 
10029 


10071 


20529 


10061 


13°62 


20635 
13630 


10C25 


19731 


AZL=A7L#+FLOAT(L)*B8(M,9)* CONS 
AZLZBZL+FLOAT (L)*B 0M, 10) *CONS 
GO TO 10030 

AZN=.90 

BZN=.0 

N2zIT +t 

BO 20929 K=1,N2 
CKPECOS(FL GAT (K-1)*H(N,1143) 
IF (K-JC) 10071,10071, 20629 
AZN=AZN-CON3*CMP*FLOAT (1) *80K, 9) FOKP 
BZN=BZN=-CON3S*CMP FLOAT (I) #90K,16) *CKP 
CONTINUE 

AZt=.0 

BZL=.6 

IF (N2-JC)109061,10030,10030 
NI=N2+4 

N0 20639 M=NLeJC 

L=M-14 
CKP=COS (FLOAT (LI FB IN, L123) 
AZL=EAZL-CONS*CMPFFLOAT (L) 490M, 9) FOKP 
C7L=BZL-CON3*CMPFFLOAT (L)*8(%,16) °CKP 
GO T) 193039 

AZN=.0 

RZN=.6 
CONL=3.1415927/SIN(B(N,11)) 
N2=Ie¢1 

0O 2cJ28 K=1,N2 
CKP=COS (FLOAT (K=1)*B(N,24)) 
IF (K-JC)10972,10072, 20028 
AZN=AZN#CONL*SMP*B(K,9)*CKP 
RZN=BZNtCON1L*SMP*B(K,10) *CKP 
CONTINUE 

AZL=.] 

BZL=.0 

IF (N2*+ JC) 10662 516030,1603C 
Ni=N2+1 

90 20338 M=N1,JC 

L=M-1 
SKPSSINC FLOAT (CL) *B(N,11)) 
AZL=AZLECONLFCMP*B(M, SP *SKP 
AZLECZLECONLFCMP*B(M,1C) *SKP 
AZCIgND=ZAZNGGZL 
BHCI»NI=BZN43ZL 

CONTINUE 

NO 10934 I=1,5C 

DO 100341 J=1, JC 


COT LI C1 e301 Sd) FCBIT 59/90, Bd -16) FAT 9G) 
CULy JPSFLOATCSD* CAZ(SeT) #8115) *BHC 5, I0) 


CALL SIMEQ(JC,C,CC) 
00 40335 I=1,5C 
B8(I,12)=.6 
8(7T,13)=.6 
B(I,14)=.6 

00 16935 J=1,J0 


ACT gL 20= BCL yl 2) FF LOAT CSP FCOC Se LD FAZC So TV/S3CT GDF (LoS (10-3 6153))) 
BC Lgl ZSP=HBCL oi SPFFLOATCS) FCCC Ue LP FBC SS TI/SIC Te GDF (Les (Le ~9 (15399) 
10335 BCTeiGI=CO(U,LI*FSINCFLOAT (J) *8 0754110) /ROT 94) +B (114) 
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<r eenmennee - 


86130416) 2.0 


20001 00 16636 I=#i.JC 


BCL 1502 (BUT eo 1G) * S(T eGdF (BIT, GIB (Te BI -AIT 423) FB CT, 49 dd MH FBS, 


12) 


B(Is16)=Bl1e159*B U1, o) 


BUTe lL 7V=Z (BCT SI/B CT eB FB UT Lh) BCI eG) ® (9015494901512) SBC Odd) Fh, 


1*8(9,2) 


BIT=ATAN(BCT. 8)) 
BTI=ATAN(B(I,5)) 


PUT 91822 20%3.14159279B1T 14) SCOS(BTIIZ(1./B(1 589 -B(T 9139) 
B19 19) =660 31% (B42) BCL 9 39581322) SZ2edFCSINIATTIS (BUT G9*9079 2) 


1*COS(BTI<-BTTI)IFF2 


it 


13°38 


13359 VCLy25) 2S TMPUN(B0193) 9801915) 5 JC) FB (692904903, 2) FFSP TL LULEI27ER (7 92 
19**2 /(166*B(992)) 
B(Ts26)=S IMPUN (8 (193) 5801516) 5 JCD FR08 9 294903 ,2) FF2*R(7 42) F354 016.% 
1 8(5,29*B8(9,2)) 


ae ee wae ~ 7 wm 


IF (3-19959e6 

IF €3-3C) 10.959 

R(T,20)=.0 

A(T,21)=.0 

GO TO 11 

CONTINUE 

BU Le 20) = (10-8 C19 7) F861 ,69 7905, 18) FB0T 2 59) FE(T 15) 
RIT ,22) 2 (20 tbls 7)*8 (156) /801 5189/8 61559) *B(1517) 
CONTINUE 
RC1e22)=39(1,20)*8 (1,4) 
CONTINUE 
PPL=ESIMPUN (3 (153298 (2516) .I5C) 
PP2=SIMPUN (8 (193) 5801515) oJC) 
PP3=SIMPUN(B (153358 (1517) 5J5C) 
PPG=PPL/PP3 

PPS=PP2/'PPZ 
SPH=SIMPUN(B (193) 9811522) JC) 
PP7=SIMPUN (8B (1,3)58(1,20),5C) 
PPB=SIMPUN (8 (1,3)e2(1,21),JC) 
PPG=PPE/PP8B 

PP1C=PP7/PP8 

NO 120039 [=1,NC 

JCIT=SC+Lel 

00 1255940 L=t, JC 
XO=B(L 23) -B(Ts 3) 

IF (XO) 860,650,871 

XO=C.0 

GO TO 861 

CONTINUE 

XO=B(L,3)~-BUL, 3) 

CONTINUE 

B2(L)=XO* BIL, 20) 
BI(LI=XO/B(Ls 3) *B IL, 21) 

IF (UCI -2)10061.,.10041,10059 
8(I,25)=.6 

R(T,26)=.0 

B(I.27)=.0 

P(I.28)=.0 

GO TO 10039 


BII=ATAN(B(1e5)) 
SBI=SIN(BTI) 


7° 
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16039 


236 


792 


TI 
734 


70¢ 


7c1 
52 
10049 


109861 
20841 


CBI=cOS(BTT) 
Bl1e¢272B (1,25) *CBI+B(l,26)*SBI 
8(1,26)=B8(1 425) *SBI-B(1,26)"C8TI 

CONTINUE 

00 206 T=190 

B(T,12)=B0I 42 7801512) 
B(T,13)=B(I 4) *B0T,13) 

BIT, 1G0=B il ge) * BC Teh) 
IF(JEE.GT.0) GO TO 10081 

615(161)=PP1C 

DO 16649 I=1.JC 

TF (B(T,6)) 702,702,703 

CC1=C.C¢ 

GO TO 704% 

CC1=B8(7518)/3(1,6) 

CONTINUE 

CLiCI)=CC1 

CC2=1.546*CC1 

CC 3=. 0679*CC1 

IF(CCL) 701,700,701 

cc4&=C.C 

#(I,36)=CC4 

GO TO 52 

CC4=B(I,7)9/CC1 

B(1,38)=CC4 

CCONS (1) =CC1 

CCTWO(I)=CC2 

CCTHR(I)=CC3 


SUBROUTINE SI MEQ(JC,CyCC) 
CCFOR(T) =CC4& : 


CONTINUE C <--- NEW VERSION AY JACK DISKIN <= X7445¢ 
RETUZN 
END DIMENSION BA (12) 90(125912) @CC(12912) 
MPO=JC 
JCI=ICH1 


CCC ICT,» 1)=0. 
O00 80 T=1,J2 
C(IeJCII=-CC(l 1? 
00 78 J=2,J/0T 
78 COCCI, JG. 
COCLe 1 FeCUL gp Itt d/CCL 54) 
8} CClIeI+i)=1. 
00 69 K=2,J/C 
00 86 T=1_MPD 
BA(I)=0. 
00 86 J=1,JCI 
66 RACID= BACT) + CCK, J)*CC II, 
MPO=MPD-1 
00 #9 Jzi, JCI 
Z=CC(19J)/84 (1) 
00 89 I21,_,MP0 
89 COI, JP=-BACT #1) 4Z ¢ COCI+1, J) 
00 99 J=z1.JCI 
99 CC(5,1P=CCIL J) 
RETURN 
END 
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SUPROUTINE OISCOT (XA, ZA, TABX, TABY, TABZ, NC pNYgN7 ANS) 
DIMENSION TASX (1), TABY (1) » TASZ (1) »NPX (37) gNOV (37), YY (37) 


C MERGE OF DISCOT,OISSERsAND LAGRAN USING BUILT IN CONSTRAINTS 


15 


2 


25 
3) 


99 


NLz=2 

IFONCT.EQe-hh) NL=Z 

ID=2*NL-2 

NUPP=NY-NL 

00 15 IIZNL,NUPP 

NLOC=II 

IF( TARX(II).GE.XAd GO TN 20 
CONTINUE 

NUPPZNUPP-RNL #2 

GO T6 99 

NUPP=NL OCN #1 

NU=NUPP FID 

00 25 JJ=ENUPP NU 

NOIS=JJ 

IF(@ TABX(JJ) -EQ. TABX(JU4L) 9 GO TO 3n 
CONTINUE 

GO TO 99 = 
NUPP=NOIS-I0 

IF (CTABXC(NOIS) eLTeXA) NUPPENDIS#! 
SUM=C.C 

NN=NUPP 

N=NN4¢I0 

00 3 I=NN»N 

PROC=TARY(T) 

00 2 J=NN,N 

A=TABX(T)-TABX(S) 

IF (A-£Q.0.) GO TO 2 

B=(XA- TABX(J) IZA 

PROO=PROD°8 
CONTINUE 

SUM=SUM#PROD 

ANS=SUM 
RE TURN 
ENO 


129 ) 


come, 
+ 


18 


10¢9 
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SUBROUTINE STRESS (AZZ, AREA, X5AR, VYBAR,AYEXI,AVYEVYOSEMXOSEMYOSEMT BoM 
10B »STRMAX,RAK, RAKOP) 

OIMENSTON CHORO(11),THICKNS(11) CAMBER (11), PITCH(11) -SKEWR E11) 
1sXI(11),0TCL11) 

COMMON/C WEIGHT /XI «CHORD, THICKNS » CAMBE ReP ITCH eSKEWR SN ITAMZZ2,0°N 
1 eRAKE sg FI ,PP7, PPB, PPO, PP11 SEWAKE pVSeRPS gSIGMAGEAR, RT P,PNO 
2 »FWODIAM,AFTOIAM,HUBLEN»FDBO°CE ,ADGORE »DISO°FL 

COMMON /UNITS/ SI,UI,UO 

OIMENSTION RAK(11) 

DIMENSION AZZ(11,38) 

OIMENSION XE(20) 

DIMENSION HAC20) ¢HAL(20) gPHI(29) gPHI2 (29) oXU1(29) 571020) ,01020),0P 
LHI (20) eSPHI (2094X%4 (20) ,AE(20),9€620),PE(29? 

OIMENSTION &(1359 ,8013),C(13) 61013) o£ C130 F013) oe GCLSDHILS) 0113), 
XP(13) 40013) C13) SCL 30 9 TELS? UCL ID GV C13) WEL TD , XEL3) eV (1L 3947613) 
DIMENSION R1(7416)5 S1(7_16) 

OIMENSION CENTST (7), CENTSO (7) 

DIMENSION FMX (7) 4 TX€7) pF MMXC 133 5 YTXO13) SKE (135) SXUCLC) 

NIMENSION VOL (7) » CENTS (7) 4 AL(79,A207) 4 X2R8°O(7) CO NTS2(7),82(13), 
x F5(13) 5P2(13),02 613) pAACLG) 589410) eCENTG2(7) »CENTMS(7), 
x TSKEWL (7), TSKEW2(7) g ASKEWL (7) ,ASKEW2 (7) 

OIMENSION V2(13) ,02013),£2(13) 

OIMENSION ALPHTA (7) 

OIMENSION XMT (10) oXL (10) 9 X4C10) »XT040796ST4(16),STLT(10) 

NIMENSION ARE A(7) ,XB8AR(7) 4 YBAR(7) ,AVEXO(7) -AYEYO(7) gEMKO(7) gTMYO(7 
1¥¢EMT93(7) gEMO9(7) gSTRMAX (7) 

DATE AF SAT ATX, ATF ATFFE pAFXX, ATTT,GRAV/1 2 40 79556720995 034152580 32 
2 9 eS 321 5 eb1733 9012714532014 / 


* ¥ * SIMPLE BEAM APPROXIMATION INCLUDING 
* © © BENDING, CENTRIFUGAL AND TORSIONAL FORCES. 
PI=3.1415926536 

ATXK=. 202084 

NN=1 
FOCPMAT (8F 946) 

ZZ=AZZ (942) 
VS=A7Z (7,2) 

DI AM=AZZ(3,2) 

DI A=CIAM 

VEL=E0.0*AZ2(5 42) 

ISEC=9 

AZZ(10,20)=0.0 

AZZ(10,21)=6T (11) 

827013 .22)=0.6 

90 100¢ I=1,10 

OCID=AZZ(1,26) 

TCIDEAZZ(1I, 21) 

EC TY=AZZ(1,22) 

C(I) =A7Z (1,23) *OTA*12.0 

SKEWCIP=ECCIIS2 60-AZZ(1924) *OIA%1 2.0/2.0 
KUCID=AZZ(1I 949) 

00 3t J=297 

AZZ(5 925) SAZZ (Jet 225) 

00 1091 I=1,7 

K=14-T 

J=K/13 

UCT eo ¥ (1-2) +.05F SF (3-1) 


UK ei-U(T) 

THCIDSAZZ(I 425) *D0 JAS12.8 
FMX(T)=.0679*DOIASAZZ(1,18) 
YSAR(ID=.0679%D0 TAP AZZ (1418) 

FMX(1)=0.0 
ACIDSCCIDSTX CT) FAT 
XCIVECCLPE ATX SAT 
YOLD)D =FMXCIIO ATE SAT/2. 

GOLISCCLIIPFTR CLI) (FE MKC I) SFQFATEF eT XC LPS SSFATIT/ 3 VHATIT SFY CID VT) 
HCTISCCIVFFSETKCIPSAT XX = ACTOFKCLTIEXIT) 
1001 CONTINUE 
C CALCULATE THE VALUE OF F1 FROM INFUT VALUES. 
26 FLA1. 99059 (DIL AMS 2.03 FF 3FVSFF2EPTFG.0/77Z 
FF U=F 1 
YFC(PD0.GT.0.) GO TO $2 
no 215 Y=1,10 
245 PCID = TALYSPLFKu(T) 
S2 NN=NN41 
CO OD GO OF CC OF 60 68 OF OF OS 6008 OD 06 06 OF O48 ON OS 08 OE OS 06 OG OE BORE COMPUTE O08 COOP OE OS OO OS OO O8 OD OO 66 OS GEES OF OF O46 04 68 000060 0000 00 046800 
CALCULATIONS FOR CONSTANTS USED IN OETERMINATION OF TOPQUS AND THRUST 
CALCULATIONS OF BENDING MOMENTS FROM THRUST AND TORNUE. 
N0 360 ITS=192 
FLsFFL 
OANL=OIAM*O.5%12.0 
IF (15-2)55,56456 
SS 00 210 T=1,1C 
PE CIY=PCT) 
RECTPSDCIDFCL HE CIDFTCID) 
BCTV SO CIP CECI VET CID) 
BE (I) =B(T) 
TUL SP CIS (OT FXUCIDS 
PHICT)=ATAN(TCTI)) 
CPHKI(LI=COS (PHICI DY 
SPHICI)=SINC(PHICID) 
HA(TI=(SKEM(ID 2 * (CPHI CT) 7 (RADNLFXUCTID)) 
YULCIV =XUCTOE COS CHA(T)) 
RAL(TDEC(IN/2. 
213 IFCSKEW(CID EQeHALCIDD KULCTIP=HKUCT) 
GO TO £8 
S56 00 57 J=1,.10 
P(I)=PECT) 
S57 301) =BE 07) 
58 F1i=F1i/EC.0 
00 69 I=1,7 
T1=I 
IF (15~-2)62,63, 63 
62 XO=XuU1L(I) 
GO TO 64 
63 xO=XUCT) 
64 I13=0 
00 68 I2=I1,.10 
T3=1I3+1 
IF (15-2) 65, 66, 66 
65 KG (I3)=(XUL(L 2) -X0) 
XE (13) =XKXU(T2) 
GO TO &€7 
66 X4(I3) © (XUCI2) -x0) 


oaqa 


Se Rae =, 1 - : * ~ Fe 


te etn 


oir 4 


AQAA 


XE (13) =XU(T2) 

67 TLCI3SV=XGC013) SAE CT 2) 

68 OLCI3=KG(I3)*90T2) 
TCT) =SIMPUN(XE , T1513) 
OCTI=SIMPUNIXE ,Q1,-13) 
TCT CL) FF FL 

69 QACID=Q(1) FFFL 


LOOP WHICH APPROXIMATES STRESS UF TO TOPSION R= SULTING 
FROM SKFW 
XTC) = LIFT FORCE » XMTC(I) = MOMENT ONE TD Lir. 


PO 111 IT=1,7 
IF (15-2) 820,830, 83C 
82> X¥MT(I)=0.69 
¥K=1.99C5*® (DI AMS 209) F¥2.FVSFF2,FPI/S 62.0877) 
KA=C(I9/2.0 
YB=TX(T)/200 
00 222 J=I49 
XL (JP SASS (SKEW (D9 ~245*C(9)) 
XL (UY =XL CL) ABS (SKEW (ID 6S *CCI)) 
XTCSPAAE CIDE LEXKS(COSCPHIC SIP CE CID) 
XM (JP=XT CUP EXEL CS) 
XMTC IT) =XMTCID XMS) 
2c2 CONTINUE 
STMCIDEXMT (IDF 2. CS(PIFXASKIFF2, ) 
STL TCIVSEXMT (IT) 42.07 (PIFXSFXAF42.0) 
GO TO A&C 
A390 STM(I)=0.00 
STLTCT Hc etl 


BGN CONTINUE 


ao 


411 CONTINUE 

559 FOPMAT (1LHO 929K oe LHK yg LOX pGHTAUM SL 2K gSHTAULS 19%, 7HM SUB T) 
S03 FOPMAT (14920X 44F24%.6) 

600 FORMAT (1H1,50X,32HSMEAPING STRESSES JUE TI TORSION) 


LOOP WHICH GALCJLATES VOL. OF SrCTIONS. 
VOLTOT=0.0 
JO 241 [T=19E 
VOL CTISACTIS (XU CT 41 P= XUCIDY*ILEM/2886 
241 VOLTOT=VOLTOT +VOL (1) 
VOL (7)2A(7)* (12 XU(7)) SO TAM/5756 
243 VOLTOT=VOLTOT+VOL (7) 
LOOP WHICH CALCJLATES CENTRIFUGAL FORCE AND STR7SS. 
TF (15-2) 264492664246 
P44 NO 245 T2196 
265 ALITY=XULCLI + C (XUL (LOL De XUL(TI D206) 
ALOT) S¥UL (7) COXUL CLO) =XUL (79972098 
GO TO 248 
246 CO 247 I=t6 
247 ALLIV=XUCT) #C CXUCT #1) -XUTID/200) 
AL(ZPEXUCT DEC CXU C10) KU7I DZ 209) 
LOGP TO TRANSFER CONSTANTS FOS DETERMINING X2°AQ, 
248 00 236 17197 
X2BAR(T) = 0.0 
236 AQUI) = ALT) * VOLED) 


C LOOP TO CALCULATE RADIAL CENTROIO ( X2BAR ). 
00 251 J=1,7 
XZAARCTID =C (A201L)94A202) A203) A264) OA2059 442 (65 4A2(7)) 7 VOLTOT 3 
x * (DIAM/2.0) 
A2(I) = 0.0 
C UNCORPECTED FORCE ANDO STRESS FOR OUTPUT OF ANSWERS WITHOUT THE EFFECT 
C RAKE AND SKEW TAKEN INTO CONSIDERATION. 
264 CENTS(I) = DENSH. OFPIFF24VELF* 2°V OL TOT* X28AQ(TI / (3600. 0*GRAV) 
CENTST(I) = CENTG(T) 7 ACT) 
251 VOLTOT = VOLTOT = VOLCT) 
C LOOKING AT THE EFFECTS OF RAKE ANO SKEW IN THE PROPFLLER. 
00 2¢3 T=1-16 
AA(T)=PIFXUCT) 
263 BBCII) = SQRTCAACT)F#2¢P(T) F¥2) 
00 267 I[=157 
TSKEWLOIT)D = (COI) 7200 - SKEW(TDY * AACTI/SB(T) 
KK = 1 
146 IF CX2BARC ID =XUCKK) FOIAM/2.0) 14945149 ,151 
151 KK> KK¢i 
IF (KK- 103146, 149,149 
149 TSKEH2¢T) = (COKKD/200 © SKEWCKK)) © AACKK)/B9 (KK) 
ALPHIACID=ATAN(TSKEW2CI) /S(X2BAR(I) * 12.0) ) 
CENTG2(IT) = SENTG(I)*COS CALPHIACT)) 
CENTMS(T) = CENTG2CTIFCTSKEW2(T) ~- TSKEWL(I)) 
CENTS2¢I) = CENTG2CT) / ACT) 
ASKEWL CLPSCTSKEWL (IDFR CI) ZAACT)) @ RAK(T) 
ASKEWZ (TIS CTSKEN2CIDFPCID/SAACIDD & RAK (KK) 
77 CONTINUE 
267 CENTMOC(I) = CENTG2CT) * ( ASKEW2CT) = ASKEWL IT) ) 
NO 281 IT=197 
JOLT) = CAT CID FCENTMOCT)) FAACTO + CQCTDOCENTYS(I)) FECT) /BBCT) 
E CIDE CIT CID OCCENTMO(ID) FP CLIN CAC I) -CENTMS(T)) *AN(TI ISBT) 
D2 CIV=(CTCIPFAACTI4QCIDFP (ID) 7 BRIT) 
281 E2CIVS(TCID FP CIP-QCUIDFAACID) 7 BBCI) 
C PROGRAM CONTINUES. 
N0 350 I=1,7 
K=1 
S(i)=2 
Z (KY=U(7) 
F2=FMX( ID +4962" TX(T) MVCT) 
00 370 L=1K 
BILI=ACCCCLIVFZC LI“ XCIID FE CIV SH CTD CF 240 (TI /GCI) tCENTS2 (1) 
RACILPZSCCCOLIFZCLIOKCI)DFESCIVISHC ID OF 2402 (1) /G CT) tCENTSTCT) 
v2(L) = ABS(B2(L)) 
300 VC(LI=ABS (BILD) 
F3=V(1) 
Fu = V2C1) 
F(TP=B(1) 
FS(I) = Bei) 
00 326 L=1,K 
FCTIEV CL) 
32d FSCLP=V2(L) 
F5S(I? = 0.0 
34D PLTDRKCIDFECTISH CTP (COLTS (LIMV CIPD FOCTISGCIDFCENTS2 (1) 
P2CTI=z-KCTIFE2CLI SHC TIM (COLD FS CLIP MOV OT) ANZ CLI/G CT) CG ENTST CT) 
O2CIIStO CTD K CLD FEQSCTIS HCI} — CoV CIPD 80207) /60T) +CENTST (TL) 
350 OCITPSCCEIVMKXCIDPFECTISHCTE mC HV CID DFOCT ISG CID ECENTS2 (FT) 
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a wae AN Ne Sag AE TE: Pe ene = 
= é~ , 
ey neg 
ae ty 7 - 


g anal n 


00 106 I#ir7 
AREACTI=ACT) 
XBARCII=KCT) 
AYEXOCT) =G6(T) 
EMXOCT)=O0¢1) 
AYEYOCT) =HUT) 
FMYO(TI=E(T) 
STRHAX(IP=AZZCOI 511) 
EMTACLIATOL) 
10) FMOB(1)=QUT) 
IF (15-2) 35153520352 
351 CONTINUE 
CALL PRNSTR CF gPeOeS TMeSTLT 9427) 
C a2e0 00000008 0000 000d es cnse 08 00 00 C000 0000 00 09 1808 OF 2090 PRINT ouTyPurT 90 6000 0000 00 0000 C000 nn ee en ee eee erence nr ee e000 00 8! 
Ge TQ 360 
352 NUMMY=DUMHY 
NN = NN#t 
360 NNENNGL 
RETURN 


( END 


rm 
i i 13 


j SUBROUTINE PRNSTR (XK VV ,729S15S52,A27) 
OIMENSTON AZZ(11 538) 
DIMENSION XXC4 CP oe V¥(100522010) eS2010),S20LCE 


CALCULATION OF PRINCIPLE STRESSES 
OUE TO TORSION ANO BENOING. 


aaNan 


NOIMENSTON XI2¢10) ,xXxI3¢10) 
00 333 K=1,7 
XT20K) =-S10K) S10) 
XT3(K) =-S27K)FS20K) 
333 CONTINUE 
XXX=0e1 
10 6446 L=4,7 
XXX=X¥XxX¢0.1 
00 555 M=1,3 
IF (M-2972422, 33 

72 XI1=KXX (0) 
OO=(ABS(KI19) **2.0 
X0=00-46,*XI2(L) 
CC=(ABS (KD) DE .5 
GO TO 4& 

22 xTi=vy tt) 
GO TO 66 

33 xTi1=ZZ(L) 

65 AC=CARSIXI1)) F#2.,0 
X¥0=DI-4.*XI3(L) 
CC=(ASS (XO) DF F,5 

44 SIGMAL=(XI1#CC0)72.0 
SIGMA2=(XT1-50172.0 
AZZ (Lo MP =EKKX 
AZZ(LeM+1C}=SIGMAL 

555 AZZ(LeM+20)25 IGMA2 
444 CONTINUE 

70) FORMAT (1He 20K 9F1262,6Xe2E 2066) 

85) FOPMAT (1HO,33X,1HX912X »GHSIGMAL »10Xg6HSIGNO?) 

860 FORMAT (1H0,2X%,99HSTRESSES AT £A4CH X STATION AOS GIVEN IN THE FOLLO 
XWING OPODER*® MIOCHORD, LEANING FCGE, TRPATLING EGS. ) 
PETURN 
ENO 


TRE, 


FUNCTION SIMPUN(X,¥,N) 
C FORTRAN IV FUNCTION FOR SIMPSONS RULE INTEGRATION 
Cc ARBITRARY NO. AND LENGTH INTERVALS K.MEALS NSRDC CODE 842 10-5-67 
OIMENSION X(2),Y (2) 
IF(N-2) 7,5 5e% 
5 S= (¥ (1) €¥ 029) © (X02) 2X01) 972. 
GO TO 6 
& M=Ne4 
8 IF (M-2) 9610,11 
11 M=M-2 
GO TO 8 
Q S=OKC20-KO199 760 FIV ALIS (Zee UX CAV XLII OX CID KILID DEV C204 (3, (X02) 
LeXCLPISOXCZIIN NK (CAVPVMNV CSD FCOUX CAD KCL EFA) SEX CID OMCLIDE CK CIV X02) 
2)))) 
L=3 
GO TO 12 
15 S=5. 
L=2 
12 M=NelL 
00 1 K=LoM,2 
S IF (A9S (X (K-10 KX 010) GE ABS (X (KD X(199) GO TO 3 
IF CARS (X (KELP X10) GTeABS(X(KID—-XCL9)09 GO TO 1L 
3 WRITE (692) Ke XCK-1) eV Koel eXCKISVCKI GN 
2 FORMATC® NON MONOTONE XC(SIMPUN) ® 93% 9 2X 91P 251204 55Xq 2E120 HS Xe IG) ous 
7 S=0. 1 
GO TO 6 3 
£1 StS¢ CX (K OLD eK (KL dD S6 FOV (KLE CZ em UX KELP MX IK ML PIS EX CKD RX (KoL DDD + 
LUV OKD FL oH OXCK ELE MX (KLIS IX CK DK (Ke LDV ECK (KD RX CK MLDS UX (K +L DR KCK) 
LIVECVCKELIF(2 Cem OX (KI MK (KL IDS CK (KEL) OX KIDD DD 
6 SIMPUN=S 
RETUPN 
ENO 
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tiny open tae stots now 


ah 


SURROUTINE WEIGHT (JC eSIGMA7,HUR,PMHC, WEIGHT 3 WO IGHTH,RAK) 


fyroarh 4 


WETGHT COMPUTES THE WEIGHT AND CENTES IF GRAVITY. THE VALUCS F 
CHORD, THICKNS, CAMBER, PITCH AND SKEWS SIMS FROM GMAIN. OIA, 
D:N, RAKE AND PT ARE SET IN STRISS. OTHEP VALUES ASE COMPHITEN 

MAKING CEPTAIN ASSUMPTIONS, 

THe ARRAY RAK USED FOR RAXE COME FOO SMAIN 


AANANAQN 


DIMENSION P(13) ,BTC129,RAKC11) 

COMMON /UNITS/ SI,UI,U0 

COMMON/CWE IGHT/X, CHORD» THICKNS, CAMBEO,PTTCH, SK= WO DT AM, Z7,0SN, CAKE 
LsPITeCTSyCPSe=PyPC eWAKE »VS,RPS,SIGMA,LEAR 0 37 9 PeP99 
2 sFWODIAMs AF TOTAM,HUBLEN»FOBOPE ,ADBORE ,OLSREFL 

DIMENSTON CHORDCLILI»THICKNS(11)-CAMBSER (41) ,PITTOH(119,SKEWP (11) 
1X11) 

NIMENSTION OISTHF (36),A(38) 

OIMENSION R(9) »PMT(9) 

DATA CNSTNT Leo CNSTNT2eCNSTNT3/. 360550 RC 71 ge 2 38/ 

OLTA THRL,THB2, THB3,THSS 46H CYLINSEXNITCAL 64 DEST e6HGNATSIS 

DATA UFE,UFS 7 SH(LBF) ,S5H(N) / 

DATA VET QUE>,USL /&HCIND HIF T) » GHOIMD = / 


****VALUES COMPUTED AND DATA OUTPUTEFF* 

THE HUB OIAMETER IS ASSUMED TO BE THE ITAMETES To THE FIOST RAD 
R&ETIO TO 8E CONSIOERED AND THE HU3 ASSUMEN TO BE CYLINDSICAL. 
THE HUB LENGTH IS ASSUMED TO EQUAL THE H'Y® JIAMET=R AND THE OTS 
THE REFERENCE LINE FROM THE BUS FACE TS TAKEN AS HALF THE HUB L 
INPUT NATA AND ASSUMED OATA WRITTEN OUT, 


ANAAAND 


CGRIX oY pHP=H¥ (X#XFZEXFVHZEV EY) J (XEXEXFVEVEYI SG, 
VOL (XV HD EH® (X*#XOX*#VFVFYI/3, 

PRES (X) VD EKEK EXE KEXEKEK EV EX EXE HEVEXEY BY SV EY SY SYSY 
ATH=THO1 

BTH=THA2 
IFC HUE.NE«0.) ATH=THB3 
IF (FUP.NE.O.) BTH=THAG 

HUBOTAM=X(1)* 0 TAM 
IF(HUS.NEeO.) GO TO 50 

DISREFL=HUBDIAM/2. 

HURL £ N=HUBOIAM 

CENGRVH=HUSLEN/2. 

GO TO 27% 

Si FWORAD=FHDDIAM/2.0 

AF TPAQ=AFTOIAM/2.0 

HUBRAD=HUBDIOM/2.0 

FRRORE=F DBORE/2.6 

ARROXE=ADBOR:/2. 


C #***HEIGHT CALCULATION**** 
273 90 1¢ I=t9JC 
40 AC(I)=CHOROC(ID* THICKNS(T) 

C We IGHT OF THE BLADES 
ASAL=SIMPUN(X, Ay JC) 
WEIGHT C=CNSINT S* CIT AM#DEN*Z7*RSAL 

( c WEIGHT OF THE HUB 
IF (HUB.EQ.9.) GO TO 20C 
: 201 HRROREZ=ARBORE +0 ISREFL* (FPS0°E-AP3ORE) /HUILEN 


138 


aa 


Cc 


200 


WSi2PI°OEN*VOL CAFTRAO,HUBRAD,OISREFL) 

MBizP IT°DEN®VOL (CARBORE ,HRBORE ,DISREFL? 
WFRisWS i -WaL 

IF (HRBORE- EQ. ARBORE) HRBOREZARBORE*. 00095901 
AFMEWS1*CGRCAFTRAD, HUBRADDISREFL) “MB 19CGR(ARBOPE ,HRBCRE »DISPEFL) 
OS=HUBLEN-OISREFL 
WS2=PI*OENFVOL CHUSRAD oF WORAD.OS) 
WA2=PIT*DENFVOL (CHRBORE»FRBORE,OS) 

WFR2ZWS 2-HB2 

F WH=O TSREFLE WER? CWS 2°CGR (HUBQAD »FWORAD,NS) 
LF (FRBORE EQ. HRBORE) FRBOREZHRBORE*. 90090901 
FWH=FWM-WB2*CGRCHRBOREs FRBORE.OS) 

WE IGHTHEWFRIOHFR2 

CENGRVH= (AFM *F WM) /SWEIGHTH 
GO T3 202 
WE IGHTH=°T *HUBOTAM*®*2*HU BLENF NENG, 


202 CONTINUE 


20 


250 
251 


WEIGHT OF THE PROPELLER 
WE IGHTP=WE TIGHT BCWEIGHTH 


*eOOCENTER OF GRAVITY CALCULATION#*##*% 
00 2c I=1s9C 
OISTHF (LY =CNST NT 2*CAMBER(TIFCOS(PITCHCI)) ¢CNSTNT 3 *CHORD(I) 


L*SINCPITCH(I)) #0 ISREFL 


THE EFFECT OF RAKE AND SKEW ARE ADDED TO THE NAYTSTANCE OF THE CE 
GRAVITY FROM THE HUB FACE FOR EACH SECTION. 
00 3C I=1eJC 
OISTHF (LI) ZOISTHF CI) “SKEWRITDFXC ID S2.FOLTAM*TAN(PITCH( I) )-PAKITI/12. 
ACT) =CHORD(T)® THICKNS (I) *OLSTHE (CT) 
AS A2Z=STHPUN(K, A, JC) 
DE NGPV3=BSA2/BSAL 
CENGFVE=DISREFL-CENGPYS 
CZNTER OF GRAVITY CONSIDERING RAKE ANN SKEW 
CFNGPV1= (WEIGHTB *CENGRVS +WEI GHTHFCENGRVH) /WEIGHTE 
CE NGPYFZNISRIFL-CENGRVI 


IF (HUB.EQ.9.) GO TO 250 
CP=PIT*DENF14.,4 
PMA=DISREFLYPMFR (HUBRAD, AFTPAD) 
PMF=D0S*PMFR(FRORAD,HUBRAD) 
PMB=HUBLEN*FPMFR (FRBORE sARRORE ) 
PHHC= CP*®(PMA+PHF-PMB) 

GQ TO 251 
PMHC=WE IGHTH* HUBRAD* #24972 

CONTINUE 

F*FPRESULTS OUTPUTF FF 
SLF=.3048 
UB=WE IGHTB 
UP=WE IGHTP 
CFL=CENGRVF/OIAM 
COL=CENGRVB/D IAM 
BA=HURLEN/OIAM 
BB=FRODIAM/DI AM 

BCZAFTOIAM/ID IAM 
BOZDISREFL/OI AH 
BEZHUBOIAM/0DI AN 


seo A 


BF=FORORE/OIA™ 
BG=ADRORE/OIAM 
SW=UFE 
IF(UO.NE.SI) GO TO 531 
SWEUFS 
U2=UB*4 448222 
IFT UP %4 2.448222 
S31 POINT tu4_SW, UB,ATH,BTH,SW,US, CFLe oil ly 
IF(HUP.EQ.0.) GO TO 55 
PRINT 106, BA,B82,B80,80,85,9F .B5 
GO TO 53 
$= CONTINUE 
PRINT 110, BE, BA, BOD 
53 CONTINUS 


MINIMUM EXPANDED ARFA PATIO CALCULATTONS 
ASS=V3/ (RPS*O0T AM) 
AJA=WAKe *AJS 
AKT=PI *CTS*AIS**2/8. 
AKQ=TPS*FAISFF3/1E. 
EARMIN=(2.646 e6*ZZ)"AKT/S(SIGMA7® (AJAFF2H(, 7FOTIFRF 254.15 
PRINT 105¢EA2MIN,AJS,AJA,AKT -AKN, oC 


124 FOPMET( 9 //25X_*WEIGHT OF BLANES*® gAS y*=¥ oF 15 64//20X 4 SWEIGHT OF P 


1ROP (BLADES +*2A6_y* HUB)* »AS,%=* 9F15e4//20Xe*CHNTER OF GRA 
2VITY OF PROP REFERENCED FROM MIDCHORND OF 2799T SECTION (= FWD, + AF 
3T) /0=* 9° 9.6//20X%,*CENTE®S OF GRAVITY JF PLADES PSFERONCED FROM 
& MIOZTHORN OF ROOT SECTION (+ FWO, + AFTI/1=* F925) 

105 FOPMETI(/20X,*KELLERS MINIMUM CAR=*,E10.% 
1//22X_*SPEED COEFF V/(ND) JS=F oF ice //2°X,* ADVANCE COFFF Ve 


21-WTT)/ (NO) JA=*,E16.4//20X,* DESIGN THONST COFFF KT=¥*,F 
310.64//20X_,*TIRQUE COEFF KQ=* E10 .4// 
G 2CX,*PROPULSIVE EFFICIENCY “TAN=*,£13.4) 

113 FOPMLT(/2CX%,*HUB DIMENSTONS/O*, BX 4* HQ OTAM =*¥F9,4/47X%, HUB L 
LeNGTH =* gFOLG/47X y*MIDCHORD OF POOT SECTION TO AFT =ND OF HU3 =*,F 
29.4) 

10> FOPMET (6F 8.4) 

107 FOFMLT( //20X,*WEIGHT OF SLADES* AS g*%y F15eb//20X%,*W TIGHT OF P 
1P0F (8LANES + TAPERFO HUBD* .A5,%=*,F 1504 4/20X%,*CENTER OF GRA 
2VITY OF PROP REFERENCED FROM MTOCHORD OF 2997 SFCTION (- FHM, + AF 
31d /=* s9FOG/S/2CKy*CENTER OF GRAVITY DF PLANTS SOCFERZ NCE) FROM 
& MIQCHORD OF POCT SECTION (- FWO, # AFT)/3=* oF9.4) 

195 FOOMAT(/20X_,*HUB NIMENSTIONS/S NAF sll X,*LINGTH=* gFALU/47KX g*FUD OTA 
LM=*%9F9.4/47K,* AFT DI AM=* 9 F9,4/G7X,*MIDCHNAN NF PNOT SECTION TO AFT 
2 ENN OF HUR=*#,F9¢G/G7K,*HUD NIAM AT MIODCHIR) OF SOOT SECTION=*,F9, 
3U/U7X_,*F HD DIAM OF BORE=* gFO,G/U7XK%,FAFT ITAY OF BORTH* F964) 

KETURN 
ce NN 
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